Isospin Physics in Heavy-ion Collisions at Intermediate Energies 



Bao-An Li^, Che Ming Kc^ and Wolfgang Baueifp 

"Department of Physics and Cyclotron Institute 
Texas A&M University, College Station, TX 77843-3366, USA 

''National Superconducting Cyclotron Laboratory 
and Department of Physics and Astronomy 
Michigan State University, East Lansing, Ml 48824-1321, USA 



Abstract 

In nuclear collisions induced by stable or radioactive neutron-rich nuclei a tran- 
sient state of nuclear matter with an appreciable isospin asymmetry as well as thermal 
and compressional excitation can be created. This offers the possibility to study the 
properties of nuclear matter in the region between symmetric nuclear matter and pure 
neutron matter. In this review, we discuss recent theoretical studies of the equation 
of state of isospin-asymmetric nuclear matter and its relations to the properties of 
neutron stars and radioactive nuclei. Chemical and mechanical instabilities as well as 
the liquid-gas phase transition in asymmetric nuclear matter are investigated. The 
in-medium nucleon-nucleon cross sections at different isospin states are reviewed as 
they affect significantly the dynamics of heavy ion collisions induced by radioactive 
beams. We then discuss an isospin-dependent transport model, which includes dif- 
ferent mean-field potentials and cross sections for the proton and neutron, and its 
application to these reactions. Furthermore, we review the comparisons between the- 
oretical predictions and available experimental data. In particular, we discuss the 
study of nuclear stopping in terms of isospin equilibration, the dependence of nuclear 
collective flow and balance energy on the isospin-dependent nuclear equation of state 
and cross sections, the isospin dependence of total nuclear reaction cross sections, 
and the role of isospin in preequilibrium nucleon emissions and subthreshold pion 
production. 
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Chapter 1 
Introduction 



Recent advances in experiments using radioactive beams with large neutron or proton 
excess have made it possible to create not only nuclei at the limits of stability, but 
also a transient state of nuclear matter with appreciable isospin asymmetry, thermal 
excitation, and compression. Extensive reviews have been written on the properties 
of these nuclei, and they can be found in Refs. 0, |^, |, ^ ^, |^, In this article 
we will review instead the role of the isospin degree of freedom in heavy-ion collisions 
at intermediate energies. The main goal of such studies is to probe the properties 
of nuclear matter in the region between symmetric nuclear matter and pure neutron 
matter. This information is important in understanding the explosion mechanism 
of supernova and the cooling rate of neutron stars. For example, the prompt shock 
invoked to understand the explosion mechanism of a type ll supernova requires a 
relatively soft equation of state ( EOS) [^J, which can be understood in terms of the 
dependence of the nuclear compressibility on isospin. In the model for prompt explo- 
]TU|, the electron-capture reaction drives the star in the latest stage of collapse to 



sion 



an equilibrium state where the proton concentration is about 1/3, which, according 
to microscopic many-body calculations, reduces the compressibility by about 30% 
compared to that for symmetric nuclear matter. Moreover, the magnitude of proton 
concentration at (3 equilibrium in a neutron star is almost entirely determined by the 
isospin-dependent part of the nuclear EOS, i.e., the symmetry energy. The proton 
fraction affects not only the stiffness of the EOS but also the cooling mechanisms of 
neutron stars |Tl|, |T2| and the possibility of kaon condensation (e~ K~Ve) in dense 



stellar matter [jT3[. If the proton concentration is larger than a critical value of about 
15%, the direct URCA process {n ^ p + e~ + u^i p + e~ ^ n + u^) becomes possi- 
ble, and would then enhance the emission of neutrinos, making it a more important 



process in the cooling of a neutron star |ri|. 

Usually, the properties of asymmetric nuclear matter are inferred from that of 
symmetric and pure neutron matter by using a parabolic approximation [|l^, and 
this seems to be supported by studies based on detailed many-body calculations. Al- 
though different assumptions have been introduced in these theoretical approaches, 
they all predict similar results for the asymmetric nuclear matter. In particular, the 
density and compressibihty at saturation are predicted to decrease as the nuclear 
matter becomes more neutron-rich. Also, the instability that leads to the liquid-gas 
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phase separation in asymmetric nuclear matter is found to arise from the chemical 
instability due to fluctuations in isospin-asymmetry instead from the mechanical in- 
stability as a result of fluctuations in baryon density as in symmetric nuclear matter 
|T^. Furthermore, the nature of liquid-gas phase transition in asymmetric nuclear 



matter is predicted to be different from that in symmetric nuclear matter. For ex- 
ample, Miiller and Serot have recently shown that the liquid-gas phase transition is 
second order rather than first order as in symmetric nuclear matter JT^|. 

To relate the theoretical predictions to the experimental observations in heavy 
ion collisions with radioactive beams, the isospin-dependent nuclear transport model, 
which uses different mean-field potentials and cross sections for the proton and neu- 
tron, has been found to be very successful. Based on this model, it has been shown 
that many observed phenomena are mainly determined by the isospin-dependent nu- 
clear equation of state and in-medium nucleon-nucleon cross sections. Other models 
such as the percolation and lattice gas models have also been found useful in under- 
standing some of the phenomena. We shall review the applications of these models 
in studying several topics of isospin physics. These include the equilibration of the 
isospin degree of freedom as a probe of nuclear stopping, the isospin dependence of 
the nuclear collective flow and balance energy, preequilibrium nucleon emissions, and 
subthreshold pion production as well as the determination of the radii of neutron-rich 
nuclei from the total nuclear reaction cross sections, An overview and outlook will be 
given at the end. 
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Chapter 2 



Isospin dependence of the nuclear 
equation of state 



Since the early work of Brueckner et al. [H| and Siemens |T7| on isospin asymmetric 
nuclear matter in the late 60's, there have been many studies on this subject based on 
different two-body and three-body forces or interaction Lagrangians. These include 
the non-relativistic Brueckner approach |jT8|, |1^, |20|, relativistic Brueckner approach 
|2l| , ^ |23| , p4| , variational many-body approach relativistic mean-field 

theory [[15], |23|, 0, 0, relativistic and non-relativistic Hartree-Fock |^ 

34, ^ or Thomas- Fermi approximation [^6[ 37 1, and the chiral sigma model [38 



Instead of discussing details of these models and their predictions, we review the 
most important, common features of the EOS of asymmetric nuclear matter predicted 
by these models. Also, we shall point out the most obvious, qualitative difference 
among the model predictions. For example, at densities greater than the normal 
nuclear saturation density, it may become important to include the relativistic effects 
39[| , which seem already to be needed in obtaining the correct nuclear matter binding 
energy and saturation density. In addition, causahty is violated at high densities 
in many non-relativistic models. We shall therefore also distinguish the predictions 
between relativistic and non-relativistic approaches. 



Fig. and Fig. ^]2| show three typical predictions for the EOS of asymmetric 
nuclear matter from the non-relativistic Skyrme-Hartree-Fock (she) model using the 
parameter set Slll the relativistic mean-field (rme) model using the parameter 
set TMl [^] and the relativistic Brueckner-Hartree-Fock approach with or without 
the momentum-dependent self-energies 1^^. The isospin asymmetry is indicated for 
each curve by the neutron to proton ratio p„/ Pp in Fig. and the relative neutron 
excess 
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P 



[2.1] 



in Fig. |2.2| . A common prediction from these studies is that the asymmetric nuclear 
matter is less stiff and bound at saturation. The minimum in the equation of state, 
i.e., the energy per nucleon versus density, disappears before the pure neutron matter 
limit is reached, and the compressibility at saturation thus decreases as nuclear matter 
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Figure 2.1: The equation of state of asymmetric nuclear matter from the Skyrme- 
Hartree-Fock (left) and relativistic mean-field (right) model calculations. The solid 
curves correspond to neutron-to-proton ratios, Pn/Pp, of (from top to bottom) 0, 0.2, 
0.4, 0.6, 0.8, and 1. Results taken from Ref. [||. 



becomes more neutron-rich. Also, the saturation density is generally reduced with 
increasing neutron/proton ratio. Although all these models give the correct saturation 
properties for symmetric nuclear matter, their predictions for the EOS of asymmetric 
nuclear matter, such as the saturation density, are quantitatively different. In the SHF 
model the saturation density depends weakly on the isospin asymmetry, while in the 
RMF model the dependence is much stronger. This difference, as we shall see in section 



2.4, results in significant differences in the density distribution and the thickness of 



neutron skins of radioactive nuclei. More detailed comparisons among results from 
the relativistic mean-field theory, relativistic Brueckner Hartree-Fock approach, and 
Skyrme-Hartree-Fock approach can be found in Refs. El| . 



2.1 Isospin dependence of the nuclear matter com- 
pressibility, binding energy, and density at sat- 
uration 

From the EOS of asymmetric nuclear matter, i.e, the energy per nucleon e(p, 6) = E/A, 
the compressibility can be calculated from 

/f(p..) = 9^. (2^2) 

5 
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Figure 2.2: Same as Fig. 2.1 from the relativistic Brueckner-Hartree-Fock calcula- 
tions. The solid and dashed curves correspond, respectively, to the case with and 
without momentum dependent self-energies. Taken from Ref. p^ . 



where the pressure P is given by 
P(p,5) = p2^e(p,<5). 



(2.3) 



Here, p and 5 are the baryon density and the relative neutron excess, respectively. In 
most models the compressibility K decreases as the matter becomes more neutron- 
rich. But the rate of decrease varies widely among models, in particular at high 
densities. 

Fig. |2.3| shows the results from a non-relativistic Hartree-Fock calculation 
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The upper and lower windows are obtained using the density-dependent Reid and 
Paris nucleon-nucleon interactions, respectively. One finds from these studies that 
neutrons tend to make the nuclear EOS stiffer. Also, at densities above about 1.5po 
the predicted density dependence of nuclear compressibility depends sensitively on the 
isospin independent part of the nucleon-nucleon interaction This is because the 
strength of isospin-dependent interaction is generally much weaker than the isospin- 
independent one. The saturation density at each relative neutron excess 5 is indicated 



by the solid dot in Fig. |2.3| and is seen to shift towards a lower density as the nu- 
clear matter becomes more neutron-rich. A similar decrease is seen for the nuclear 
compressibility at saturation density as a result of its strong density dependence. 

A non-relativistic study of the isospin dependence of the nuclear compressibility 
and density at saturation has been carried out by Kolehmainen et al. using the 
Skyrme-type interaction in an extended Thomas-Fermi approach. Their prediction 
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Figure 2.3: The nuclear compressibility K obtained from the non-relativistic Hartree- 
Fock calculations using the density- dependent Reid (left window) and Paris (right 
window) nucleon-nucleon interactions, respectively. Results taken from Ref. p5| . 




Figure 2.4: Nuclear compressibility as a function of relative neutron excess, obtained 
from the extended Thomas-Fermi calculations using three different Skyrme interac- 
tions. Results taken from Ref. |]36 |. 
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Table 2.1: Parameters Kq and a for the isospin dependence of the compressibihty at 
saturation. 



Force Paris SKM* Si' Slll AVl4+uvil uvl4+uvil FP 



Ko{0) 185 216.6 370.3 355.3 



209 



202 



238 



2.027 1.988 1.272 1.275 



2.196 



2.049 



1.01 



Ko(\) 143.3 168.7 318.0 305.0 



158 



156 



212 



on the isospin dependence of the nuclear compressibility at saturation Kq{6) is shown 
in Fig. p.4| , and their results can be parameterized as a quadratic function of 6 



Ko{6)=Ko{0){l-a6'), (2.4) 

For small S this parameterization works well and has been confirmed by other stud- 
ies. For example, Bombaci and Lombardo pOf obtained a similar result using a 
non-relativistic Brueckner-Bethe-Goldstone approach with the Paris interaction. We 



summarize in Table 2A the parameters Kq{0) in units of MeV and a extracted from 



these studies and those from fitting the EOS calculated by Wiringa et al. 



2% using the variational many-body (vmb) approach with the Argonne and Urbana 



two-body potentials AVl4 and UVl4 augmented by the Urbana three-body potential 



uvil, as well as that from Friedman and Pandharipande (fp) are also listed. As 
one expects, both -ft'o(O) and a extracted from different models or the same model 
but different nuclear interactions vary appreciably. 

We also show in Table 2.1 the compressibility -ft'o(|) at 5 = 1/3, which is about 
20-30% below Kq{0) and is of special interest in nuclear astrophysics. According 
to the prompt explosion model of supernova 0, a neutron- rich system is formed 
during the collapse of massive stars as a result of neutrino trapping, and the isospin 
asymmetry 5 of a star stays almost at a constant value of 1/3. In this study, Eq. 
( [2. 41 ) with the parameter set of SKM* was used, and it was found that the softening of 
nuclear EOS due to the isospin asymmetry is responsible for generating a successful, 
prompt explosion. To our best knowledge, this is probably the only astrophysical 
evidence that supports the predicted isospin dependence of nuclear compressibility 
at saturation. In laboratory experiments, a 30% reduction in compressibility seems 
too small to have a measurable effect even on nuclear collective flow which is believed 
to be the most sensitive probe to the compressibility. However, the corresponding 
changes in the mean-field potentials for neutrons and protons are opposite in sign, and 
this makes certain observables, such as the neutron to proton ratio of preequilibrium 
nucleons in heavy-ion collisions, sensitive to the isospin dependence of the EOS. We 



shall discuss this in more detail in section K2 
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Figure 2.5: Relative compressibility at saturation obtained from the Dirac-Hartree- 
Fock calculations using isoscalar mesons only (dashed line), both isoscalar and isovec- 
tor (solid line), and the non-relativistic Hartree-Fock calculations using the Skyrme 
Slll interaction (dot-dashed line). Results taken from Ref. 
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There are also many calculations of the nuclear compressibility using a relativistic 
approach. Although it is not easy to compare different model calculations as they 
usually involve many parameters and assumptions, the comparison made in Ref. p3| 



is very instructive. Shown in Fig. |2.5| is a comparison of the isospin dependence of 
the ratio Ko{5) / Ko{0) calculated in the framework of a Dirac-Hartree-Fock approach 
using only isoscalar {a, uj) mesons (curve b) and both isoscalar and isovector (tt, p) 
mesons (curve e) as well as from the non-relativistic Hartree-Fock approach using the 
Skyrme sill interaction (curve Slll). The general trend of decreasing compressibility 
with increasing isospin asymmetry 6 is common to all three models. The rate of 
decrease is, however, sensitive to the details in the model. 

Other properties at saturation, such as the binding energy, pressure and density, 
also show dependence on the isospin asymmetry S of nuclear matter. Shown in Fig. 
are the binding energy and density at saturation as functions of 5^ P3[. As for 



the compressibility at saturation, the saturation density can be parameterized as 

Peg(5) = Po(0)(l-M'), (2.5) 



In Table |2.2| we summarize the values of the parameters b and po(0) extracted from 
several models. Results from the Dirac-Hartree-Fock approach are denoted by dhf (b) 
and DHF (e) for using isoscalar mesons only and both isoscalar and isovector mesons, 
respectively. From tables |2.2| and ^TT] one can clearly see a large variation of the 
parameters describing the isospin dependence of the compressibility and saturation 
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Figure 2.6: Same as Fig. 2.5 for the energy (left) and density (right) at saturation. 
Resuhs taken from Ref. [P^. 

density. Obviously, experimental data are needed to test these models and thus to 
narrow down the isospin-dependence of the nuclear EOS. In this respect, collisions 
of neutron-rich or radioactive nuclei are promising tools. We shall come back to this 
point later. 



Table 2.2: Parameters po and b in the isospin dependence of saturation density. 



Force Paris SKM* Si' Slll dhf (b) dhf (e) 



Po(0) 0.289 0.1603 0.1553 0.1453 0.1484 0.1484 
b 1.115 0.634 0.286 0.084 0.65 0.9 
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2.2 Empirical parabolic law and the nuclear sym- 
metry energy 

For asymmetric nuclear matter at densities away the saturation density, various the- 
oretical studies, e.g., Refs. |TT], |2^, ^ have shown that the energy 
per nucleon can be well approximated by 

e(p, 6) = Tf{p, 6) + Voip) + S^V^ip), (2.6) 

where 

is the Fermi-gas kinetic energy, Vo(p) and ^(p) are the isospin-independent and - 
dependent potential energies. Higher-order terms in 6 are negligible. For example, the 
magnitude of the 5^ term at po has been estimated to be less than 1 MeV fl^, |26| . 
Eq. (p.6|) can be further approximated by 

e(p,5) = e(p,0) + esy^(p)5^ (2.8) 

where e(p, 0) is the EOS of symmetric nuclear matter and 

esyM = ^^^^1^=0 = Itf{p,0) + V,{p) (2.9) 

is the bulk symmetry energy. At normal nuclear matter density, the two terms in the 
above equation have about similar magnitude. Eq. ( p.8| ) is known as the empirical 
parabolic law for the EOS of asymmetric nuclear matter and is considered to be valid 
only at small isospin asymmetries. However, many non-relativistic and relativistic 



calculations have shown that it is actually valid up to 5 = 1. In Fig. |2.7| and Fig. 



two examples from the non-relativistic Brueckner-Bethe-Goldstone approach |^ and 



the relativistic-Brueckner-Hartree-Fock approach [Q, respectively, are shown for the 
total binding energy as a function of 5^ at several densities p. In both cases, the fit 
using the parabolic law shown by solid lines is indeed found to be valid in the whole 
range of 6. 

Using the empirical parabolic law, one can easily extract the symmetry energy 
esym(p) from microscopic calculations. According to Eq. (|2.8| ) the bulk symmetry 
energy esym(p) can be evaluated from the two extreme cases of both pure neutron 
matter and symmetric nuclear matter via 



^sym 



(p)=e(p,l)-e(p,0). (2.10) 



Furthermore, the symmetry energy can be expanded around the normal nuclear mat- 
ter density po, i.e.. 



11 




Figure 2.7: Total binding energy per nucleon as a function of 5^ for densities (from 
top to bottom) of 0.4, 0.3, 0.17, and 0.076 fm~^, obtained from the non-relativistic 
Brueckner-Bethe-Goldstone approach. Sohd hues are the fits using the parabohc law 
to that from numerical calculations. Taken from Ref. [pO . 




Figure 2.8: Same as Fig. 2.7 obtained from the relativistic Brueckner-Hartree-Fock 
approach. The lines correspond to values of the density of (from top to bottom) 0.18, 
0.14, 0.1, 0.06, and 0.02 fm'^. Results taken from Ref. 
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Table 2.3: The slope and curvature of symmetry energy in units of MeV at normal 
nuclear density. 

Force Paris SKM* Si' Slll dhf (b) dhf (e) 



68.8 45.78 35.34 9.91 132 138 



Ksy^ 37.56 -155.9 -259.1 -393.7 466 276 



In the above, L and Ksym are the slope and curvature of the symmetry energy at 
normal density, i.e., 

L ^ 3po^|p=po, (2.12) 

Ksyra = 9Po g2p \p=Po- (2-13) 

With this expansion, the predicted symmetry energy at normal nuclear matter den- 
sity can be compared with that determined from the mass formula, which gives a 



value of esym(Po) in the range of 27-36 MeV Most models are tuned to give an 
esym(po) within this range. For example, in the non-relativistic Hartree-Fock theory 
||46| , 0, the predicted value of esym(po) is between 27 and 38 MeV depending on 
the nuclear interactions used in the calculation. However, the relativistic mean-field 
(rmf ) theory p8|, |3^, |5^, ^ usually gives higher values of esym(po) in the 



range of 35-42 MeV. What distinguishes these models around normal nuclear matter 
density are thus the slope L and curvature -ft'sym- This is clearly illustrated in Fig. 
|2^ where results from several non-relativistic approaches |19|, ^ are compared. 
Although all predict a esyin(/9o) in the range of 27-30 MeV (in agreement with that 
from the empirical mass formula), it is seen that the slope and curvature at po are 
very different. As shown in Table p.3| , -ft'sym ranges from positive 466 MeV to negative 
400 MeV. The most recent data from giant monopole resonances indicate that -ft"sym 
is between —566 ± 1350 MeV to 34 ± 159 MeV depending on the mass region 
of nuclei and the number of parameters used in parametrizing the compressibility of 
finite nuclei. The Kgym data thus still cannot distinguish these model calculations. 

In radioactive nuclei as well as neutron stars, the symmetry energy at densities 
away from the normal nuclear matter saturation density is relevant. Although the 
Bethe-Goldstone theory predicts a p^^^ dependence for the symmetric energy, the 



density dependence is linear in the rmf theory |28|, ^ Q . Results from most other 
calculations cannot be simply described. In the most sophisticated calculations by 
Wiringa et al. using the variational many-body (vmb) theory |2^, the density de- 



pendence of esym(p) has been calculated using either the Argonne two-body potential 
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Figure 2.9: A comparison of the density dependence of the symmetry energy predicted 
by several non-relativistic many-body calculations: Paris (solid line), AVl4 (dotted 
line), uf14 (dot-dashed line), and Paris+3BF (dashed line). Results taken from ref. 

m MM- 




Figure 2.10: Symmetry energy predicted by variational many-body calculations: 
AV14-FUVII (solid line), UVl4-Fuvil (dot-dashed), and uv14-|-tni (dashed). Results 
taken from Ref. |2^. 
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AVl4 or Urbana uvl4 together with either the three-body potential uvil or TNI, and 
the results differ appreciably as shown in Fig. \2.10\ With three-body interactions the 
symmetry energy is seen to vanish at high densities. Therefore, pure neutron matter 
was predicted by the VMB theory as the ground state of dense nuclear matter. This 
prediction has important consequences on the structure and magnetic properties of 
neutron stars as stressed recently by Kutschera et al. 
predict a nonzero symmetry energy at high densities. 



54]. However, other models 
Among all models, the pre- 
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Figure 2.11: Contribution from nuclear interactions to the symmetry energy in the 
relativistic mean field (rmf) theory (solid line) and the variational many-body (vmb) 
theory (dashed line). Results taken from Ref. ||5^ 



dictions from the RMF theory and the VMB theory using the UVl4-(-TNI interaction 
differ the most. Since the kinetic contribution to the symmetry energy is about the 
same in all model calculations, the differences among model predictions are mainly 
from the V2{p) term in Eq. ( |2.9| ). As shown in Fig. |2.11| the V2(p) term is always 



repulsive and increases linearly with density in the rmf theory, while in the vmb 
theory it changes from repulsion to attraction as the density increases. The variation 
of V2(p) with density in the vmb theory can be understood in terms of the behavior 



of nuclear interactions in dense nuclear matter |^5|]. At high densities the short-range 
repulsion dominates and is greater for a nucleon pair in an isospin singlet (T = 0) 
than in an isospin triplet (T = 1) state. Pure neutron matter is therefore more sta- 
ble. At moderate densities the strong attractive isospin singlet tensor potential and 
correlations keep the isospin singlet pairs more bound, and symmetric nuclear matter 
is thus more stable than pure neutron matter. These features do not exist in the rmf 
theory where the asymmetry energy is due to the rho meson exchange, which leads 
to a repulsive V2(p) at all densities. As a result, these two approaches predict very 
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different structures and chemical compositions for neutron stars 

It is clear from the widely diverging model predictions that despite great theoreti- 
cal efforts our knowledge of the EOS of asymmetric nuclear matter is rather limited. In 
particular, the behaviour of the symmetry energy at high densities is most uncertain 
among all properties of dense nuclear matter 



2.3 Neutron stars and the EOS of asymmetric nu- 
clear matter 



The mechanism for supernova explosion and the properties of neutron stars have 
been subjects of much interest; for recent reviews, see e.g., Refs. |60 |. 
Here we only concentrate on the effects of the isospin-dependence of the EOS, in 
particular the density dependence of symmetry energy, on the properties of neutron 
stars. Various studies have indicated that the symmetry energy affects mainly the 
chemical composition of neutron stars [jll], |3^, |6^. Other properties, such 



as the cooling rate, lepton profiles and neutrino flux, which depend on the chemical 
composition, will therefore also be affected. 




Figure 2.12: Equilibrium fraction of protons in neutron stars predicted by the rela- 
tivistic mean- field (rmf, solid line) and variational many-body (VMB, dashed line) 
theories. Taken from Ref. [Q. 

In the most simple picture, a neutron star is composed of neutrons, protons and 
electrons with a proton fraction of 

x = l{l-6). (2.14) 
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The condition for /^-equilibrium in terms of the chemical potentials of electrons (fie), 
neutrons (/i^) and protons (/ip) is 

Aie = /^n - /ip = — = 4esym(p)(l " 2x). (2.15) 

For relativistic degenerate electrons of density pe = Pp = xp, charge neutrality requires 

fie = {ml + plf'^ ^ hc{3n^pxY/\ (2.16) 

which together with Eqs. ( p.l5| ) and ( |2.8| ) determine an equilibrium proton fraction 
X given by 

hc{37c^pxY/^ = 4esym(p)(l - 2x). (2.17) 

The equilibrium proton fraction x is therefore determined solely by the nuclear sym- 
metry energy, esym(p). At high densities such that p,e > both electrons and muons 
are present at /3-equilibrium. The inclusion of muons mainly alters the value of the 
equilibrium proton fraction, x, but not its density dependence. Therefore, the differ- 
ence in X predicted by using different symmetry energies is about the same with or 
without including muons |n|, As discussed in the previous section, the RMF and 



the VMB theory with the AVl4-(-TNI interaction differ most in their predictions on the 
symmetry energy at high densities. Their predicted equilibrium proton fractions at 
high densities are therefore also very different as shown in Fig. p.l2| . We see that the 
RMF theory predicts a linear increase, while in the VMB theory the proton fraction 
in neutron stars gradually decreases as the density increases. The disappearance of 
protons in neutron stars is a common feature of the VMB theory, although the critical 
density at which this occurs depends on the interaction used in the calculation. As 
a result, there exists an instability in neutron stars with respect to the separation 
of protons and neutrons at low densities. This has thus led to the suggestion that 
polarons, which are localized protons surrounded by neutron bubbles, can be formed 
in neutron stars |Q . In Fig. p.l3| we compare the equilibrium proton fractions (lower 



window) corresponding to different symmetry energies (upper window) 

esyn,(p) = (2^/3 - 1)^EUu'/' - F{u)] + e,y^{pMM)F{u), (2.18) 

with F{u) having one of the following three forms 
2^2 

Fi{u) = —-, (2.19) 
1 + u 

F2{u) = u, 
Fsiu) = u^l\ 

where u = p/ p^i^ the reduced baryon density and i?^ is the Fermi energy. These forms 
of the symmetry energy resemble closely the three typical results from microscopic 
many-body calculations discussed in the previous section. Although in all three cases 
X increases with density, the differences among them are still appreciable. 
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Figure 2.13: Left window: symmetry energy as a function of density. Right window: 
corresponding equilibrium fraction of protons. Solid lines: F2{u), dashed lines: F^lu), 
dash-dotted lines: FUu). Results taken from Ref. p!T||. 



The variation in symmetry energy, especially at high densities, results in a corre- 
sponding variation in the proton fraction x. This has significant effects on the cooling 
of neutron stars and the kaon condensation in dense stellar matter. In the so-called 
standard model for neutron stars, their cooling is mainly due to the modified URCA 
process 

{n,p) + p + e~ ^ {n,p) + n + Ue, {n,p) + n {n,p) + p + e~ + i7e- (2.20) 
The direct URCA process 

p + e~ + i7e, p + e~ ^ n + Ue 



n 



(2.21) 



is usually forbidden by energy-momentum conservation. However, it has been shown 
by Lattimer et al. |]TT| that if the proton fraction is higher than a critical value of 
about 1/9, the direct URCA process can occur. This would then enhance the emission 



of neutrinos; so the neutron star cooling rate is increased significantly. From Fig. |2.13 
it is seen that the symmetry energy determines entirely whether the proton fraction 
can exceeds the critical value and at what density this happens. 

Another important effect of the symmetry potential is on the kaon condensation 
in dense stellar matter. Since electrons with chemical potentials higher than the kaon 
mass will trigger the process e~ K^u^, the condition for this to happen is deter- 
mined by the equation /ig > m^-. It was shown in ref. [|13] that the critical density 
for kaon condensation depends sensitively on the form of the symmetry potential. 
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Figure 2.14: Equilibrium fraction of protons as a function of density obtained from 
the relativistic mean field theory. Solid curve: Qp = 5.507, dashed curve: Qp = 2.78. 
Results taken from Ref. 



The dependence of equilibrium proton fraction on the underlying nuclear interac- 
tion can be illustrated by the recent work of Sumiyoshi et al. using the relativistic 
mean-field (rmf) theory p3|, The symmetry energy is found to vary almost 



linearly with density, and its strength is related to the p meson-nucleon coupling 
constant Qp via 



esym(po) = ; : ^ + ^?;Po/2m^, (2.22) 

where M* is the nucleon effective mass, rup is the p meson mass, and kp is the Fermi 
momentum. The first and second terms are the kinetic and potential contributions 
to the symmetry energy, respectively. From the above expression and Eq. ( p.l7| ), 
one obtains the Qp dependence of the proton fraction as shown in Fig. 2.14. As Qp 



increases from 2.78 to 5.507 the proton fraction is seen to increase by about a factor 
of two. 

Eq. ( p.22| ) also shows that the effective mass of a nucleon, M*, affects the kinetic 



contribution to the symmetry energy. As nuclear matter becomes heated, which is 
expected in neutron stars, the nucleon effective mass may change. Donati et al. |^ 
have recently studied the coupling of the mean-field single-particle motion to the sur- 
face vibration within the quasiparticle random-phase approximation and found that 
M* decreases with increasing temperature for T < 2 MeV, which are the temper- 
atures relevant for the pre-supernova collapse of massive stars. The corresponding 
increase of the symmetry energy was estimated to be about 2.5 MeV/nucleon, which 
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would then increase the shock energy by about half of the observed explosion energy 
of SN1987a. The outcome of the explosion from detailed numerical simulations thus 
depends sensitively on the symmetry energy. However, more recent studies of Dean 
et al. , using the Monte Carlo shell model and taking into account pairing correla- 
tions, indicate that the symmetry energy changes by less than 0.6 MeV at the neutron 
star temperature. The temperature dependence of the nuclear symmetry energy thus 
remains controversial. 
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Figure 2.15: Mass of a neutron star M/Mq (in solar mass units) as a function of 
its central density for various proton fractions, x. Solid line: x = 0, dot-dashed: 

0- 



X 



0.15, dashed: 



X 



0.3, dotted: x = 0.45. Results taken from Ref. 



The structure or density profile of a neutron star with an isotropic mass distri- 
bution can be determined by solving the well-known Tolman-Oppenheimer-Volkov 
equation, i.e.. 



dP _ [p(r) + P(r)] [M(r) + 47ir^P{r) 
dr — 2rM{r) 



(2.23) 



where P(r) is the pressure obtained from the EOS via Eq. ( p.3|) and M(r) is the 
gravitational mass inside a radius r. In Fig. |2.15| and Fig. |2.16| the mass and radius 
of a neutron star from a relativistic Brueckner Hartree-Fock approach are shown as 
a function of central density for various proton fractions [^. A maximum mass 
of about 2.4 Mq and a radius of about 12 km is obtained, and this is rather large 
compared to the observed value of about 1.9 Mq, indicating that the underlying 
EOS is too stiff. Pion and kaon condensates as well as effects of hyperons have been 
suggested as possible mechanisms to soften the EOS |6^. Large variations of both 
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Figure 2.16: Radius of a neutron star as a function of its central density for various 
proton fractions. Solid line: a; = 0, dot-dashed: x = 0.15, dashed: x = 0.3, dotted: 
X = 0.45. Results taken from Ref. ||65[| . 



neutron star mass and density profile are observed by varying the proton fraction. 
However, the structure of a neutron star is indirectly affected by the symmetry energy 
through the equilibrium proton fraction. As shown in Refs. the neutron 

star structure changes only slightly by varying F{u) or gp in the ranges discussed 
above, and the stiffer the EOS is the less important is the symmetry energy. This 
is mainly because the symmetry energy is much smaller compared to the energy per 
nucleon in symmetric nuclear matter. Nevertheless, as we have mentioned earlier, the 
stiffness of the nuclear EOS is reduced by the neutrons. This reduction is important 
for generating a successful, prompt supernova explosion M. 



2.4 Radioactive nuclei and the EOS of asymmetric 
nuclear matter 

The radius, thickness of neutron skins, deformation, binding energy, density distribu- 
tions, and other properties of radioactive nuclei near the drip lines depend sensitively 
on the isospin-dependence of the nuclear EOS. Tanihata has recently stressed the 
possibility of extracting the EOS of asymmetric nuclear matter through the investiga- 
tion of these properties In this review we discuss a complementary approach in 
extracting the isospin-dependent EOS via reactions of neutron-rich nuclei at interme- 
diate energies. We shall first discuss the dependence of the properties of radioactive 
nuclei on the isospin-dependent EOS and then demonstrate how these properties are 
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Figure 2.17: Density distributions for the nuclei ^^epb, "^^^Fh, ^2°Sn, ^^Zr, and ^^^Ca as 
predicted by the Skyrme-Hartree-Fock (left) and the Relativistic Mean-Field (right) 



model. Results taken from Ref. 40 



In Fig. |2.1| the predictions on the EOS from the Skyrme-Hartree-Fock and the 
relativistic mean-field theories are compared. The properties of neutron rich nuclei 
have also been studied using these models. In Fig. p. 17] , density distributions predicted 



by the two models for several nuclei are shown A comparison of the results for 
^°^Pb and ^^^Pb is particularly interesting. One sees that, although the central density 
from SHF does not change from 208p]^ ^^^Pb, it decreases in rmf. This is related to 
the fact that in SHF the saturation density essentially does not change, but in rmf it 
decreases rapidly as nuclear matter becomes more neutron-rich as shown by the thick 



dashed lines in Fig. |2?1|. As a result, the rmf gives a more prominent neutron skin 
for 266pb than for ^^^Fh. 



In Fig. |2.18| , we compare the quadrupole deformation, f32, of the ground states 



for sulfur isotopes predicted by the two models |67]. The deformation for even-even 
isotopes 28-38g jg fairly similar in both SHF and rmf models, i.e., prolate {j32 > 0) 
ground states in 28-32,38g^ oblate minimum in ^^S, and spherical shape in the magic 
nucleus ^^S with N=20. However, there are significant differences in heavier isotopes, 
in particular, for ^°'^^S. They are predicted by the rmf to have prolate ground states 
with ^2 ~ 0.25; the oblate minima with (32 ~ —0.16 lie about 4 MeV above the ground 
state. In the SHF calculations, however, prolate (/?2 ~ 0.25) and oblate {P2 ~ —0.24) 
configurations are essentially degenerate. 

A recent experiment at Michigan State University by Scheit et al. has con- 
firmed the predicted large deformation of 40,42g^ ■\^^^ could not distinguish the sign 
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Figure 2.18: Quadrupole deformations of sulfur isotopes predicted by the SHF (dashed 
hne connecting the squares) and RMF URCA (sohd hne connecting the circles) models. 
Results taken from Ref. [|6^ . 
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Figure 2.19: Comparison of the quadrupole deformations of sulfur isotopes. Points 
with error bars are the experimental data, while the dashed, solid, and dotted lines 
are the predictions of the SHF, RMF and shell model, respectively. Results taken from 
Ref. |8|, ||. 
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of deformation. In this experiment the energies and B{E2;0^^ — > 2^^) values of 
the 2i states of 38,40,42 g ^gj^-g measured by Coulomb excitations. In Fig. |2.19| the 
measured deformation |/32| is compared to the predictions of the SHF, rmf and shell 
model. Both, SHF and rmf calculations slightly underpredict the deformation, while 
the shell model of Brown slightly overpredicts the deformation. Clearly, additional 
information about the sign of deformation is needed to distinguish the prediction of 
these models about the nuclear EOS . 




Figure 2.20: Binding energies of Ca isotopes. Solid points are the experimental data, 
while results from the Landau- Vlasov calculations with various values of c (c = 28: 
solid line, c = 20: dashed line, and c = 0: dot-dashed line) are shown as diamonds 



connected by lines. Results taken from Ref. [70 



Properties of neutron-rich nuclei have also been studied using the Boltzmann- 
Uehling-Uhlenbeck (buu) or Landau- Vlasov (lv) models fT^, 0. These models, as 
we shall discuss in more detail in section [4.1| , were originally developed for studying 
the reaction dynamics of heavy-ion collisions. They are based on the semi-classical 
approximation, and are thus suitable for studying the dependence of the gross prop- 
erties of neutron rich nuclei on the isospin-dependent EOS. 

To include the effects of isospin on nuclear dynamics, Farine et al. have added to 
the Zamick potential energy density, which has been widely used in standard buu/lv 



models, an asymmetric term of the form |70 



tfsym = Cp[{pn - Pp)/po]^ (2.24) 

where the coefficient c = esym(po) ~ is the the symmetry energy at normal nuclear 
matter density due to nuclear interactions. They have studied the dependence of 
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Figure 2.21: Same as Fig. 2.20 for the Fermi energies of protons and neutrons. Taken 
from Ref. 



ground state properties of neutron-rich Ca isotopes on the isospin-dependent EOS by 



varying the value of c. Fig. |2.20| shows a comparison of the calculated binding energies 
of three neutron-rich Ca isotopes with the experimental data shown by plus signs. The 
data for ^^Ca do not exist, its binding energy is obtained using the extended Thomas- 



Fermi with Strutinsky- integral model (ETFSI) of Pearson et al. As expected, 

nuclei become less bound with increasing c. Moreover, this effect is the strongest for 
the most neutron-rich nucleus ^^Ca. However, the model tends to overbind nuclei 
for c up to 28 MeV. The comparison of the Fermi energies of neutrons and protons 
for the three Ca isotopes are shown in Fig. p.21| . General features of the Fermi 



energies are seen to be correctly described by the LV model. The results from c = 
are clearly far from the data since the Coulomb energy of protons is then not 
properly balanced by the symmetry energy. As a result the density distribution of 
neutrons and protons obtained with c = is incorrect. Since the symmetry potential 
is repulsive for neutrons and attractive for protons, a stronger symmetry potential 
would deplete more neutrons from the center to the surface of a nucleus as noticed 
in Refs. |^ |7^, |75[. The demonstration by Sobotka using the BUU model and 
an isospin-dependent EOS based on the Skyrme ll effective interaction is instructive 
[71[| . The density distributions for neutrons and protons in ^^''Au obtained from this 



study are shown in Fig. p.22| for the case with (solid line) and without (dashed line) 



symmetry potential. Shown are (from top to bottom) the total nucleon density, the 
neutron density, and the proton density. In both cases, the total nucleon density 
distribution is the same. However, there are differences in the individual density 
distributions for protons and neutrons. Without the attractive symmetry potential 
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for protons, they are pushed towards the surface by the Coulomb force, and there is 
almost no neutron skin. Including the symmetry potential, a neutron skin of about 
1/3 fm is seen, and this is consistent with the droplet model prediction 
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Figure 2.22: BUU model calculations of the neutron and proton densities for ^^'^Au 
with (solid lines) and without (dashed lines) symmetry energy. Shown are (from top 
to bottom) the total nucleon density, the neutron density, and the proton density. 
The calculations are from Ref. ||71||. 



To summarize, properties of neutron-rich nuclei are sensitive to the the EOS of 
asymmetric nuclear matter at saturation density. However, to extract properties 
of isospin asymmetric nuclear matter at densities away from the saturation density 
collisions of neutron-rich nuclei will be more useful. Since semiclassical transport 
models can reproduce the gross properties of neutron-rich nuclei and their dependence 
on the isospin-dependent EOS, they thus provide a useful framework for studying 
collisions of neutron-rich nuclei and for extracting information about the isospin- 
dependence of the nuclear EOS. 



2.5 Chemical and mechanical instabilities in hot 
asymmetric nuclear matter 

Nuclear matter is not thermodynamically stable at all densities, p, temperatures, T, 
and neutron excesses, 6. The necessary and sufficient conditions for the stability of 
an isospin-asymmetric nuclear matter can be expressed by the following inequalities 
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dT 



> 0, 



p,6 



(—] > 



OS 



> 0, 



(2.25) 
(2.26) 
(2.27) 



P,T 



where E, P and yU„ are the energy per nucleon, pressure and neutron chemical po- 
tential, respectively. The first condition ensures the thermodynamic stability and 
is satisfied by any reasonable equation of state; the second condition prevents the 
mechanical instability against the growth of density fiuctuations; and the last one 
protects the diffusive stability against the development of neutronization. The last 
two conditions might be violated in certain regions in the (p, T, S) configuration space. 
Nuclear multifragmentations initiated by various mechanical instabilities, such as the 
volume instability, the surface instability of the Rayleigh kind, and the Coulomb in- 
stability, have been extensively studied during the last decade [^, |8l|, |8^; for 



a recent review, see, e.g., [Q. With the recent advance in radioactive ion beams, 
a nuclear matter with a large isospin asymmetry can be created transiently in these 
reactions. It is thus of interest to study the chemical instability of isospin-asymmetric 
nuclear matter. 
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Figure 2.23: RMF predictions of pressure as a function of density at fixed temperature 
T for various proton fraction x (left window), and at fixed x for various temperature 
T (right window). Taken from Ref. |T6[| . 
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To understand the respective roles of chemical and mechanical instabilities on the 
dynamics of nuclear matter, it is important to study the relative locations of these 
instabilities in the (p, T, 6) configuration space. Since the adiabatic spinodal is always 
enclosed by the isothermal spinodal (its), it is sufficient to compare the diffusion 
spinodal (ds) and ITS. Obviously, ;the boundaries of these instabilities depend on the 
equation of state. Using RMT theory Miiller and Serot have found that the diffusive 
spinodal encloses more of the configuration space than the isothermal spinodal |]TB|. 
Fig. |2.23| shows typical results form their calculations at a constant temperature T 
(left window) or a constant proton fraction x (right window). At currently reachable 
region of x > 0.3 in neutron- rich nuclei induced reactions, the separation between the 
isothermal and diffusive spinodals in P — p configuration space is appreciable. 

As discussed in previous sections, the equation of state of asymmetric nuclear 
matter varies significantly among models. In particular, the rmf theory predicts a 
characteristic linear density dependence of the symmetry energy. It is therefore of 
interest to examine how the conclusions based on the rmf might change when other 
equations of state for asymmetric matter are considered. In the following we discuss 
the results for the chemical and mechanical instabilities in a thermal model using 
phenomenological equations of state. The advantage of a thermal model is that at 
modestly high temperatures, most quantities relevant to the discussion of chemical 
and mechanical instabilities can be evaluated analytically. 



300 
% 200 
X 100 

0.' 



to 

-100 

a. 



> 
<u 



Oh 



-3(D0 
5 




-10 



1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 

; (5=0.0 

: F(u)=u 


1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 

(a) - 




T=4 (MeV) : 
T=8 (MeV) " 
T=12 (MeV) - 


\" " ' " ' 

L\ 

:-.\ 

\ " ^ 


/-/: 

/./ _ 

/ -/ : 
/ ./ 


A " ^ 


/ -/ 

^ - / 


1 1 1 1 1 1 1 1 1 





1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — 1 — r 

; (5=0.2 ; 


"I — 1 — 1 — 1 — 1 — 

(b) ■ 






\"- \ y -'/ '■ 


■ \ "- ^ — ^ .-'/ 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



0.0 0.2 0.4 0.6 



0.8 0.0 

P/Po 



0.2 0.4 0.6 0.8 
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ditions as functions of density 
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=4, 8 and 12 MeV for b =0.0 
1871. 
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Figure 2.25: Same as Fig. 2.24 for 6 =0.4 (left panel) and 6 = 0.6 (right panel). 
Taken from Ref. P7[] . 



The density Pg of neutrons (g = n) or protons (g = p) is 
1 

Pg = ^ k^fg{k)dk, 

TT^ JO 



where 



fq{k) = [exp{eq - Pq)/T + 1] ^ 



(2.28) 



(2.29) 



is the nucleon Fermi distribution function, and Cg is the single particle kinetic en- 
ergy. At modestly high temperatures (T > AMeV), this equation can be inverted 
analytically to obtain the chemical potential [0, |86| 



IJ^q = Vq + T 



n=l 



n 



(2.30) 



where 



27ch 
niT 



2 \ 1/2 



(2.31) 



is the thermal wavelength of the nucleon, and h'^s are inversion coefficients. In Eq. 
( p.30|) Vq is the single particle potential energy and can be parameterized as 



Vq{p, s) = a{p/po) + Kp/poT + vnp, 5). 



(2.32) 
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The parameters a, b and a are determined by the saturation properties and the 
compressibihty K of symmetric nuclear matter, i.e., 



-29.81 
23.45 



46.90- 



K + 44.73 



K 

is: + 255.78 



K- 166.32 
K + 44.73 



- 166.32 
(MeV), 



(MeV), 



211.05 



(2.33) 
(2.34) 
(2.35) 



The isospin-independent term should also contain a momentum-dependent part 
which is important for certain dynamical observables in heavy ion collisions, such as 
the collective flow (e.g., [^, 0), but is not essential for observables such as the 
neutron to proton ratio of preequilibrium nucleons Effects of the momentum- 

dependent interactions on thermal properties have been recently studied in Refs. 
p^ , p3| , and rather small effects on the global properties of nuclei are found. Here we 
neglect the momentum-dependent interaction and shall concentrate on investigating 
the effects of the isosp in- asymmetric potential V^^^. This potential is given by 



VUp,S)=dwaip,5)/dp, 



(2.36) 



where Wa{p,S) is the contribution of nuclear interactions to the symmetry energy 
density, i.e.. 



and 



Wa{p,5) = Ca ■ pF{u)5'^ 



ea = esym(piVA^) " (2^^^ - 1) l^j 



Using the three forms of F{u) discussed in section i.e., F{u) 
the corresponding symmetry potentials are, respectively. 



u 



and 



yn(p) 
asy 

yn(p) 
asy 



ynip) 
asy 



2r2 



±2eaU 6 + CaU 6 



±2eau6, 



(2.37) 

(2.38) 

u and M^/^, 

(2.39) 
(2.40) 

(2.41) 



With the chemical potentials determined by Eq. ( |2.30| ), it is then possible to 
obtain the total pressure of the nuclear matter by using the Gibbs-Duhem relation 



dP 
dp 



;i + 5)%^ + (i-5) 



dp.. 



dp dp 
The pressure can be separated into three parts, i.e., 

P = -Pkin + -Pq + -Pasy) 



^2.42) 



(2.43) 
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where Pkin is the kinetic contribution 



P 



kin 



Tp 



n=l ^ 



and Po is the contribution from the isospin-independent nuclear interaction 
1 



Pn 



Po 



(2.44) 



:2.45) 



Pasy is the contribution from the isospin-dependent part of the nuclear interaction, 
and can be written as 



asy 



asy 



2e.po(-)''5', 
Po 

eaPo(— ) , 
Po 



and 



P 



asy 



2 Po 



(2.46) 
(2.47) 

(2.48) 



for F(u) = u^, u and u^^^, respectively. 

From the chemical potential and pressure, the stability conditions can be deter- 
mined. The determination of the mechanical stability condition is straightforward. 
However, to calculate dnn/d6 at constant pressure P and neutron excess 6, the fol- 
lowing transformation is useful: 



06 



^ 06 

'OfXn 

06 



^OP ^ 

' Opn 

. 9p , 



T,5 



T,5 



'0P\ 
^06) 

'0P\ 
dp 



T,p 
T,5 



'OP' 
06 



(2.49) 
(2.50) 



T,P 



Knowing {Ofin/06)rp p and {OP/Op)rp ^ as functions of p, T, and 6, regions of 
chemical and/or mechanical instability can then be identified in the configuration 
space. First, we show in Fig. |2.24| and Fig. |2.25| these two quantities as functions 
of p at various T and 6. In this calculation K = 200 MeV and F[u) = u are used. 
For comparison results for symmetric nuclear matter are shown in the left window of 
Fig. |2.24| . In this case there is no chemical instability and the quantity {Onn/06)j, p 
increases with both density and temperature. To understand this, we note that 



ii) 

Opn \ 

06 



T,p 



^^asy + ^ + l)bn{^r [(1 + " (1 " 



(2.51) 



T,P 



'Qyn ^ 
asy 

06 



T,P 



+ T 



n+l 



(2.52) 
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and 



lim 



5^0 




0. 



Then, it is easy to show that 

) T,P 



hm. 



(2.53) 

(2.54) 
(2.55) 



which is positive as one expects for symmetric nuclear matter and increases with both 
density and temperature. 

The mechanical instability is known to occur at intermediate densities, i.e., in the 
mixed phase, between the gas and liquid phases at subcritical temperatures. This is 
clearly demonstrated by plotting {dP/dp)rp ^ as a function of p. It is interesting to 
examine this quantity in the limit of low densities. From the Gibbs-Duhem relation 
of Eq. ( |2.42| ) and the expression for pg of Eq. ( p.3(]| ), one has 




lim 



T,5 



(T + a-^) 
Po 



(2.56) 



as expected for an ideal gas. 

We note from Eq. ( |2.49| ) that the quantity {dpn/dS)rp p is singular along the 
boundary for mechanical instability, where {dP/dp)rp ^ = 0. This singularity is indi- 
cated by the spikes in Fig. |2.24| and Fig. p.25| . However, their heights at the boundary 
of the gas and mixed phases is very small because the last derivative in Eq. (|2.49|) is 
close to zero, i.e., 



limp^o ^ 




(2.57) 



Furthermore, in the ideal gas limit the nuclear matter is also chemically stable, i.e.. 



lim 



T 



'dpr, 



dS J 1 + 5 



> 0. 



(2.58) 



As the system becomes neutron-rich, the chemical instability starts to develop, 
so the boundary for mechanical instability shrinks while that for chemical instability 
expands as 6 increases. For example, in contrast to the variation of the chemical 
instability the mechanical instability at a temperature of T=10 MeV gradually disap- 
pears as 6 increases from to 0.6. Furthermore, in the mechanically unstable region 
the system is chemically stable, but at higher densities chemical instability can ap- 
pear in the mechanically stable region. Therefore, in neutron-rich nuclear matter 
the chemical instability is a more relevant quantity in determining its properties. As 
shown in Fig. ^1^, where the pressure along the chemical, (ds), and mechanical, 
(its), spinodals are plotted as functions of density, both the isothermal pressure and 
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Figure 2.26: Pressure as a function of density at fixed temperatures T =4, 8, 12 and 
16 MeV for 6 =0.0, 0.2, 0.4 and 0.6. The isothermal and diffusive spinodals are shown 
by the dashed and dotted hues, respectively. Taken from Ref. |^ . 
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Figure 2.27: Chemical (upper window) and mechanical (lower window) instability 
conditions as functions of 5 at a fixed temperature T = 10 MeV and various densities. 
Taken from Ref. ||87|1 . 
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Figure 2.28: Pressure as a function of density at a fixed temperature T = 10 MeV 
at various 5. The isothermal and diffusive spinodals are shown by the dashed and 
dotted hues, respectively. Taken from Ref. [p7| . 



the separation between the two spinodals increase with increasing 5, so not only the 
chemical but also the mechanical instability is strongly isospin-dependent. 

The isospin-dependence of mechanical and chemical instabilities at a fixed tem- 
perature are studied in more detail in Fig. |]2^ and Fig. 1^281 . It is seen that nuclear 
matter is both mechanically and chemically stable at low and high densities. At inter- 
mediate densities (e.g., p/po = 0.4, 0.5) the system is mechanically unstable for small 
5 (e.g., 5 < 0.2 and 5 < 0.4 for p/po = 0.4 and 0.5, respectively.), and chemically 
unstable for intermediate S (e.g. 0.2 < 6 < 0.4 for p/po = 0.4 and 0.4 < 6 < 0.6 
for p/po = 0.5, respectively.) The corresponding boundaries in the pressure-density 
space are shown in Fig. p.28| . Again, the diffusive spinodal is more extended than 
the isothermal spinodal. The dependence of pressure on isospin is also shown in this 
figure. Its strong dependence on 6 is due to the significant Pasy contribution to the 
total pressure. 

The above analysis thus indicates a strong isospin-dependence of both chemical 
and mechanical instabilities. It is therefore of interest to know how this dependence 
might be affected by the isospin-dependent and -independent parts of the nuclear 
equation of state. This is illustrated in Fig. p.29| and Fig. |2.30| , where the pressure is 
shown as a function of density (left panels) and 6 (right panels) along the diffusive 
spinodals (upper windows) and isothermal spinodals (lower windows) at a constant 
temperature of T = 10 MeV by using the three forms of F{u) and a compressibility 
of 380 and 200 MeV, respectively. Several interesting observations can be made 
from these two figures. First, the diffusive spinodals are always more extended in 
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Figure 2.29: Pressure as a function of density (left panel) and isospin asymmetry 
(right panel) at a fixed temperature T = 10 MeV along the boundary of diffusive 
spinodals (upper windows) and isothermal spinodals (lower windows) by using the 



three forms of F(u) and a compressibility of K = 380 MeV. Taken from Ref. ^1 
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Figure 2.30: Same as Fig. 2.29 for a compressibility of K = 200 MeV. Taken from 
Ref. mi. 
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the configuration space than the isothermal spinodals. Second, both DS and its 
depend on the form of F{u); and this dependence is stronger for the stiff equation of 
state with K = 380 MeV. Third, from the right panels of these two figures one sees 
that it is both chemically and mechanically favorable for nuclear matter to be less 
asymmetric (smaller S) in the liquid phase than in the gas phase. This result can also 
be obtained from energy considerations. Since the equation of state for asymmetric 
nuclear matter contains a esym(p)^^ term, it is therefore energetically favorable for 
asymmetric matter to separate into a liquid phase that is less asymmetric and a gas 
phase that is more asymmetric, rather than into two phases with equal asymmetry. 
Finally, there is a maximum asymmetry for both DS and its. As the density decreases 
from the liquid phase the asymmetry increases along both DS and its to a maximum 
value. The asymmetry then decreases towards the gas phase. As shown in these 
figures, the maximum asymmetry is also sensitive to both the isospin-dependent and 
-independent parts of the nuclear equation of state. 

The results discussed in this section are more than just of academic interest. Sig- 
nificant differences have been observed in the multifragmentation of isospin-symmetric 
and isospin-asymmetric nuclear matter formed in heavy ion collisions at intermediate 
energies ^ |^. BUU calculations suggest that the abundance of neutron 



rich isotopes observed in experiments is primarily the result of the coalescence of the 
clusters with neutrons from the gas phase That these neutron exist in abun- 

dance is a direct consequence of the large isospin-asymmetry in the gas phase. It is 
interesting, by the way, that the percolation model has difficulties in explaining the 
fragmentation data for neutron-rich isotopes [Q, despite its success in matching a 
large number of other oberservables in nuclear fragmentation. 

One popular explanation for nuclear multifragmentation is that the hot nuclear 
matter formed in the reaction expands adiabatically into the mechanical instability 
region where it fragments into clusters and nucleons due to the growth of density 



fluctuations ||9^. In describing this process great efforts have been devoted to incor- 
porating effects of density fluctuations into dynamical models during the last decade 
100| , |101| , |102| , p.03| , |104|| . To describe multifragmentation in asymmetric nuclear col- 



lisions using these models, it will be necessary to extend them to include the isospin 
degree of freedom and its fluctuations. This will then allow one to compare the 
strengths of both density and isospin fluctuations and to study the relative time 
scales of chemical and mechanical instabilities in asymmetric nuclear matter. Indeed 
some interesting new features about the collective motions new the n/p drip lines and 
the onset of chemical and mechanical instabilities have been found in recent studies 
using RPA||105[| and Landau dispersion relation approaches | P- 6|] . 



2.6 New features of liquid-gas phase transition in 
asymmetric nuclear matter 

The liquid-gas phase transition in asymmetric nuclear matter is not only more com- 
plex than in symmetric matter but also has new distinct features. This is because the 



37 



nature of phase transitions in a matter is strongly influenced by its dimensionality. 
The inclusion of the isospin degree of freedom in studying the phase transitions in 
nuclear matter has thus attracted much attention recently [jT6|, |5^, |78|, |107| |108 



109| , p.lO|| . These studies have shown that there are indeed new features associated 



with the liquid-gas phase transition in asymmetric nuclear matter. 




First of all, the critical temperature for the liquid-gas phase transition has been 
predicted to decrease with increasing neutron excess. This can be understood qual- 
itatively in terms of the increasing contribution from the asymmetric pressure Pasy 
Fig. |2.31| shows the critical temperature as a function of neutron-excess 6 predicted 
by the isospin lattice gas model of Kuo et al. [|110|| . This prediction is in qualita- 
tive agreement with that based on the Skyrme interaction and rmf theory ||16|| . 
Using the latter theory, Miiller and Serot have studied the sensitivity of the critical 
point to the variation of the symmetry energy esym(po) and the bulk compressibility 
K. In Fig. p.32| the binodal sections predicted by the rmf theory at T = 10 MeV 
for esym(po) =30, 35 and 40 MeV are shown by the solid, dotted and dashed curves, 
respectively. As one would expect from our discussions in the previous section, the 
minimum proton fraction at which the phase separation occurs increases only by 
10%, and the critical temperature decreases by less than 1% as the symmetry energy 
increases from esym(po) =30 to 45 MeV. On the other hand, the critical pressure in- 
creases by about 40%. However, the shape of binodal surfaces are rather similar for 
different values of symmetry energy. As we have stressed earlier, in the rmf theory the 
symmetry energy has a characteristic linear density dependence. Other models with 
different density dependence in the symmetry energy naturally may predict critical 
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Figure 2.32: Binodal sections at T = 10 MeV and esym(po) =30, 35 and 40 MeV 
shown by the sohd, dotted and dashed curves, respectively. Taken from Ref. |T^. 

points at different temperature, pressure and proton fraction. It is well-known that 
the critical point of liquid-gas phase transition is sensitive to the bulk compressibility 
K. Miiller and Serot have found that the critical temperature changes by about 13% 
when K varies from 200 to 300 MeV using the same symmetry energy esym(po) = 35 
MeV. Although this variation is larger than that due to the variation of the symmetry 
energy, it is unlikely to have a significant effect on the liquid-gas phase transition. 

An important question about the liquid-gas phase transition in asymmetric nu- 
clear matter is its order. Glendenning ||109|| as well as Miiller and Serot |TB[ have 



recently stressed that the isospin degree of freedom plays an important role in phase 
transitions. In the liquid-gas phase transition, not only the total baryon density but 
also the total charge density is conserved. For an asymmetric matter, the latter con- 
straint leads to a smoother variation of the thermodynamical quantities, such as the 
temperature, specific entropy, and heat capacity, during the phase transition. The or- 
der of phase transition in asymmetric matter is then expected to be higher than that 
in symmetric nuclear matter. In Fig. |2.33|, the variation of temperature and specific 



entropy for proton fractions of 0.5, 0.45, 0.4 and 0.3 are shown. In the symmetric case, 
the specific entropy increases by 3 units in a temperature interval of 1.5 MeV, which 
increases to about 5 MeV as the proton fraction decreases to 0.35. Thus, both the 
specific entropy and temperature change continuously during the phase transition, 
so the heat transferred to the system is used not only to convert the liquid to vapor 
but also to heat the matter. The first-order liquid-gas phase transition in symmetric 
nuclear matter therefore changes to a second-order one in asymmetric nuclear matter. 
Another important question about the liquid-gas phase transition in asymmetric 
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Figure 2.33: Specific entropy as a function of temperature at a constant pressure for 
proton fractions of 0.5, 0.45, 0.4 and 0.35, respectively. Taken from Ref. |T^. 



matter is the relative proton fraction in the two phases. Since the EOS of asymmetric 
matter contains a term Csym^^, it is energetically favorable for the system to separate 
into a neutron-rich liquid phase and a proton-rich gas phase rather than two phases 
with equal proton concentration. The same conclusion has been reached from the 
study of chemical and dynamical instabilities in the previous section. However, one 
should be cautious in comparing this conclusion directly with experimental observa- 
tions on the charge/mass ratio of heavy fragments and single particles or light clusters 
because the liquid-gas phase transition may occur in heavy residues during the later, 
expansion phase of nuclear reactions. In the latter case, the charge/mass ratio of the 
heavy-residue depends strongly on the symmetry potential during the early stage of 
the reaction. We shall discuss this in more detail in section [4.2| . 

To confirm experimentally these new features of liquid-gas phase transition in 
asymmetric matter is challenging. One possible way is to study the dependence of 
nuclear multifragmentations on the charge/mass ratio of the reaction system. 
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Chapter 3 



Isospin dependence of in-medium 
nucleon-nucleon cross sections 



In extracting information about the structure of radioactive nuclei, such as the radii 
and the distributions of the constituent neutrons and protons, the model used most is 
Glauber's. In the optical limit, the total interaction cross section is determined by a 
transmission function, which is a convolution of the nucleon-nucleon cross section and 
the density distribution of nucleons from both the target and projectile in the over- 
lapping region. Knowing the proton density distribution (which can be determined 
from other means such as electron scattering) and the nucleon-nucleon cross sections, 
the neutron density distribution can then be determined. Usually, only the isospin 
averaged free nucleon-nucleon cross section is used as input. The effects of including 
the isospin-dependence of the in-medium nucleon-nucleon cross section needs to be 
studied. Also, the isospin-dependent in-medium nucleon-nucleon cross sections are 
needed in transport model to extract the isospin-dependent EOS and to understand 
the isospin-dependent phenomena found in heavy-ion collisions. In this chapter we 
will first discuss the isospin dependence of free-space nucleon-nucleon cross sections 
which are normally used as inputs in the isospin-dependent transport models. Then 
we shall review recent theoretical studies on the in-medium nucleon-nucleon cross 
sections based on the many-body theories discussed in the previous Chapter. 



3.1 Isospin dependence of free-space NN cross sec- 
tions 

It is well-known that the scattering cross section between two nucleons depends on 



their isospin. Fig. |3.1| compares the free-space cross sections for neutron-proton and 
proton-proton or neutron-neutron scattering as functions of bombarding energy. The 
data in the energy range of 10 MeV < -Eiab < 1000 MeV can be parameterized by 
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= -70.67 - 18.18/3-^ + 25.26/?-^ + 113.85/5 (mb), (3.1) 
cjf;"^ = 13.73 - 15.04/3-1 + 8.76/3^"^ + 68.67/?^ (mb), (3.2) 
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where /? = f/c is the velocity of the projectile nucleon. It is seen that below about 
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Figure 3.1: Cross sections of neutron-proton and proton-proton scatterings as func- 



tions of bombarding energy. Taken from Ref. ||112 



500 MeV the neutron-proton cross section is about a factor of 2 to 3 higher than 
the proton-proton or neutron-neutron cross section, indicating that the T matrix for 
the total isospin 7 = channel is much stronger than that of the 1 = 1 channel. In 
nuclear medium, one expects that the strength of G matrix, which differs from the 
free space T matrix because of the Pauli effects and many-body effects, for these two 
channels are also different. How their relative strengths may change at finite density 
and temperature is a question of current interest. Since the reaction dynamics of 
heavy-ion collisions is governed by both the nuclear EOS and the in-medium nucleon- 
nucleon cross sections, both are important and have to be studied simultaneously. 
In fact, the main uncertainty in the extracted nuclear EOS from heavy-ion collisions 
is due to our poor knowledge of the nucleon-nucleon cross sections in medium. It is 
thus desirable to find experimental observables that are only sensitive to either the 
EOS or the cross sections. This is particularly important for studying the symmetry 
potential since it is quite weak compared to the nuclear EOS of symmetric nuclear 
matter. As we shall discuss in the next chapter, transport models are very useful for 
such studies. 
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3.2 Isospin dependence of in- medium NN cross 
sections 



The study of in-medium nucleon-nucleon cross sections have been mostly based on the 
Bethe-Goldstone (bg) equation fS^, |g |m|, |TT|, |TT7|, [TTB|, [Tg. Medium 

effects appear in the BG equation mainly through the Pauli blocking factor for inter- 
mediate states and the self-energies of the two nucleons in the denominator of the 
propagator. However, results from these studies differ appreciably, with some models 
predicting a decrease of nucleon-nucleon cross sections while others an increase. 

In the Dirac-Brueckner approach of Refs. |117|| , model parameters are fixed by 
fitting free-space nucleon-nucleon scattering data and deuteron properties. Nucleon- 
nucleon cross sections in nuclear medium at zero temperature have been predicted to 
decrease with increasing density. For example, at normal nuclear matter density and 
a bombarding energy of 50 MeV, both o"„p and reduced by about a factor of 



two. Results of the calculations in Ref. ||117|| can be parameterized by 
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Fig. |3.2| shows the neutron-proton and proton-proton cross sections at normal and 
twice normal nuclear matter density. One sees that the strong isospin dependence of 
nucleon-nucleon cross sections in free-space disappears gradually as the bombarding 
energy and density increase. Furthermore, at densities larger than about twice the 
normal density, the proton-proton cross section even becomes higher than the neutron- 
proton cross section. The respective effects of the Pauli blocking and selfenergy 
corrections have not been studied in Ref. | 117 |. Bohnet et al. |TT^| have also studied 
the in-medium nucleon-nucleon cross sections during the collisions of two slabs of 
nuclear matter at zero temperature. The Pauli blocking factor has been estimated 
using two Fermi spheres separated by the beam momentum, and it is found that the 
in-medium cross section generally increases with density. 



In the work by Aim et al. [ 118 , 119], the cross section in a hot nuclear matter has 
been evaluated by including also the hole-hole collisions in the Pauli blocking operator. 
Fig. ^]3| shows their results at temperatures of 10 MeV and 35 MeV and a density of 
0.5po as a function of bombarding energy. Effects of the Pauli blocking and selfenergy 
corrections are separated by comparing full calculations with those by setting the 
Pauli blocking operator equals to 1. First, it is seen that at both temperatures the 
selfenergy correction suppresses the cross section, while the Pauli blocking operator 
for intermediate states enhances the cross section. Second, at energies above about 
200 MeV predictions for different temperatures converge to values smaller than the 
free-space cross section. It is also seen that at lower temperatures the cross section 
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Figure 3.2: The neutron-proton (solid lines) and proton-proton (dashed lines) cross 
sections as functions of bombarding energy at normal (thick lines) and twice normal 
density (thin lines). The calculations are taken from Ref. ||117|| and are compared to 
the free cross sections of Fig. 3.1 (solid and open circles). 



has a strong peak above the free-space cross section. This has been interpreted as a 
precursor effect of the superfluid phase transition in nuclear matter ||118| , |119|| . 



In summary of this section, our current theoretical understanding about the in- 
medium nucleon-nucleon cross section, especially its dependence on isospin, density 
and temperature, is very limited. Most calculations are for nuclear matter only. 
Ideally, for heavy-ion collisions one should evaluate the in-medium cross sections 
selfconsistently at densities and temperatures determined by the reaction dynamics, 
and this requires much more work in the future. It is worth mentioning that recent 
studies on collective flow in heavy-ion collisions at intermediate energies have shown 
an indication of reduced in-medium nucleon-nucleon cross sections ||12CI| , |121| , |122 



An empirical in-medium nucleon-nucleon cross section ||121 



^medium /i \ „ P \ Jree /o ci\ 

(^NN =(!+" — Jf^WTV (3-6) 

Po 



with the parameter a ~ 0.2 has been found to better reproduce the flow data [|120 



122 . It is thus very promising that heavy-ion collisions may shed some light on the 



in-medium cross sections and their isospin dependence. 
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Figure 3.3: Isospin averaged nucleon-nucleon cross section as a function of bombard- 
ing energy at a density of O.Spo- Taken from Ref. ||lim . The dotted line represents the 
free cross section, the sohd hues represent the in-medium cross sections at tempera- 
tures of 10 (thin hne) and 35 MeV (thick hne). The dashed hues are the corresponding 
cases without Pauh-blocking. 
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Chapter 4 

Isospin-dependent phenomenology 
in heavy-ion collisions at 
intermediate energies 



Several new and interesting phenomena related to the isospin have been observed 
in heavy- ion collisions with radioactive beams ||123 |. To understand such phenom- 
ena many existing models developed for heavy-ion collisions, such as the isospin- 
dependent percolation model [Q, isospin-dependent lattice gas model ||1 10|| and 
isospin-dependent transport models, have been extended to include the isospin degree 
of freedom. We shall discuss in this chapter some of these phenomena and their expla- 
nations. Transport models are particularly useful for studying the isospin-dependent 
EOS and the in-medium nucleon-nucleon cross sections, and we shall therefore also 
discuss in this chapter the role of isospin in these models and their predictions. 



4.1 Isospin-dependent transport models for heavy- 
ion collisions 

A powerful framework to test various predictions on the isospin dependence of nuclear 
EOS and in-medium nucleon-nucleon cross sections is the nuclear transport theory. 
The isospin degree of freedom has been included in transport models at various levels 
of complexity J70|, |7l|, 0, 0, |T2|, |125|, |T2|, [I23, |12|, g |13|. These models were 
developed from the well-known Boltzmann-Uehling-Uhlenbeck (buu) transport the- 
ory [Ml, |13|, [13|, |l3l or the Quantum Molecular Dynamics (qmd) theory [|13|]. by 
evolving separately the phase space distribution functions for protons and neutrons. 
In these models, the initial phase space distribution of nucleons in the target and 
projectile, nucleon-nucleon cross sections, mean-field potentials, and Pauli blocking 
are all isospin-dependent. For the isospin-dependent nucleon-nucleon cross section, 
the free-space cross sections shown in Fig. have usually been used. 

For the isospin-dependent mean-field potential, equations of state based on nuclear 
many-body theories discussed in section p.2| have been used. For illustration we 
discuss here the properties of symmetry potential given in Eq. ( p.39|) . For the simplest 



46 



form of F{u), i.e., F{u) = F2{u) = u, one has 



K:?^ = ±2e,u6 = ±2e, 



Pn - Pp 

Po 



(4.1) 



This is the asymmetric part of nuclear mean-field potential used in Refs. |73| , p.26| , |129 
136| , |137|| . The asymmetric energy density of Eq. 
to the same mean-field potential at p ~ po a^s Eq. 



used in Refs. fT^, leads 
^Jl). To illustrate the magnitude 



and variation of the symmetry potential, we show in Fig. H7 
three forms of F{u) given in Eq. ( p.l9|) and esym(po) 



V;:(p)(p,(5) using the 



32 MeV. It is seen that the 




Figure 4.1: Symmetry potentials for neutrons and protons corresponding to the three 
forms of F{u) (see text). Taken from Ref. [|13(]| . 



repulsive (attractive) mean-field potential for neutrons (protons) depends sensitively 
on the form of F{u), the neutron excess 6, and the baryon density p. In collisions 
of neutron-rich nuclei at intermediate energies, both 6 and p can be appreciable in 
a large space-time volume where the isospin-dependent mean-field potentials, which 
are opposite in sign for neutrons and protons, are strong. This will affect differently 
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the reaction dynamics of neutrons and protons, leading to possible differences in their 
yields and energy spectra. 

In modeling heavy-ion collisions the nuclear mean-field potential should also in- 
clude a Coulomb term Vq for protons, i.e., 

\/"(^)(p, 5) = a{p/po) + bip/poT + + KS'Hp, 5)- (4.2) 

To extract esym(p) from the experimental data requires observables that are sen- 
sitive to the asymmetric part but not the symmetric part of the nuclear EOS. This 
is necessary because the magnitude of V^y^-* is rather small compared to the sym- 
metric part in Eq. ( [4.21) . In addition, these observables should not depend strongly 
on other factors that affect the reaction dynamics, such as the in-medium nucleon- 
nucleon cross sections. In this respect, it is interesting to mention that the possibility 
of extracting the isospin-dependent EOS from the elastic scatterings of neutron-rich 



nuclei has been examined recently by using a double-folding model p^, |138|| . Both 
the isovector and isoscalar parts of the interaction potential between two neutron-rich 
nuclei are evaluated. 
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Figure 4.2: The symmetry and isoscalar potential for the scattering of ^^Li C at 
a beam energy of 60 MeV/nucleon. Shown is the calculation utilizing the BDM3Y1- 
Paris potential. The dashed line represents the symmetry potential, and the solid 
line is the isoscalar potential. The results are taken from Ref. |^ . 



Shown in Fig. are the isovector and isoscalar potentials for the scattering of 
^^Li+^^C at a beam energy of 60 MeV/nucleon. The strength of the symmetry poten- 
tial (isovector) is much smaller than that of the isoscalar potential, and consequently 
the optical potential and thus the elastic scattering cross section is dominated by 
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the isoscalar potential as shown in Fig. O. As a result, it is difficult to extract the 
symmetry potential from elastic scattering data. 




Figure 4.3: Elastic scattering cross section for ^^Li C at a beam energy of 60 
MeV/nucleon. Left window: the solid line is the result of a full calculation; the 
dashed line is obtained without the symmetry potential; and the dash-dotted line is 
obtained using the real potential folded with unrealistic compact density distribution 
for neutrons in ^^Li. Taken from Ref. 
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These studies further indicate the importance of selecting appropriate observables 
in order to learn the isospin-dependent part of the nuclear EOS. As in charge-exchange 
reactions, which are useful for extracting the isospin-dependent part of the nucleon- 
nucleon interaction ||139|, 1140 



the ratio of preequilibrium neutrons to protons in 
heavy-ion collisions at intermediate energies have been found to be suitable for ex- 
tracting the isospin-dependent EOS [^. We shall discuss this ratio in detail in the 
next section using the isospin-dependent transport model of Refs. fT^, |127| p.28| , |129 



To describe properly reactions induced by neutron-rich nuclei, it is important to 
use as realistic as possible initial neutron and proton distributions in the transport 
model. This can be achieved by initializing neutrons and protons randomly in a sphere 
with a radius corresponding to the rms radius of a heavy nucleus, as we have discussed 
in section . Then the BUU / LV model with a proper symmetry potential can describe 
the gross properties of the nucleus after evolving it in time, despite the fact that 
the model does not contain any information on its microscopic structure. Since the 
binding energy as well as density distributions of neutrons and protons depend on 
both the strength and form of symmetry potential, i.e., one obtains different density 
profiles using different symmetry potentials, it is thus possible to learn about the 
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asymmetric EOS from studying saturation properties of neutron-rich nuclei. However, 
the information inferred from these studies is only around the saturation density. 
What is relevant for astrophysics, but poorly known, is the asymmetric EOS at other 
densities. These can only be reached in heavy-ion collisions. The effects of different 
symmetry potentials on particle emission in heavy ion collisions might be stronger 
than on equilibrium properties of neutron-rich nuclei. To observe these effects the 
same initial states for both the target and projectile have to be used in calculations 
using different symmetry potentials. In Refs. [p.29| , |130|| , this has been achieved by 
using initial nucleon distributions given by the rmf theory. 



4.2 Isospin dependence of preequilibrium nucleon 
emission 

In this section, we will first discuss the ratio of preequilibrium neutrons to protons 
{Rn/p{E-kin) = dNn/dNp) from collisions of neutron-rich nuclei at intermediate ener- 
gies by using the isospin-dependent BUU model. Then we will discuss experimental 
observations. 




r (fm) 



Figure 4.4: Density distributions of neutrons (dotted) and protons (solid) in Sn iso- 
topes calculated using the RMF theory. 
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As an example, collisions of ^^^Sn +^^^ Sn, ^^"^Sn +^^^ Sn and ^^^Sn +^^^ Sn at a 
beam energy of 40 MeV/nucleon are studied. The first two reactions have been re- 
cently investigated experimentally at nscl/msu by the MSU- Rochester- Washington- 
Wisconsin collaboration |^4|, Preequilibrium particles have been measured in 
these experiments and are now being analyzed ||141|| . The last system is included 
for discussions and comparisons. The density distributions of neutrons and protons 
in the Sn isotopes are calculated using the rmf theory | |142| | and are shown in Fig. 
4^ , and the local Thomas-Fermi approximation is then used to generate the nucleon 
momenta. To identify free nucleons, a phase-space coalescence method has been used 
at 200 fm/c after the initial contact of the two nuclei, when the quadrupole moment 
of the nucleon momentum distribution in the heavy residue is almost zero, indicating 
the approach of thermal equilibrium. A nucleon is considered as free if it is not cor- 
related with other nucleons within a spatial distance of Ar = 3 fm and a momentum 
distance of Ap = 300 MeV/c. Otherwise, it is bounded in a cluster. These results 
are not sensitive to these parameters if they are varied by less than 30% around these 
values. 

Effects of the compressibility K of symmetric nuclear matter and the in-medium 
nucleon-nucleon cross sections on the ratio i?„/p(i?kin) can be studied by neglecting 
both the Coulomb and symmetry potentials in the BUU model. In Fig. |4.5| this 
ratio is shown as a function of nucleon kinetic energy for central (upper window) 
and peripheral (lower window) collisions of ^^^Sn +^^^ Sn at a beam energy of 40 
MeV/nucleon. The total number of preequilibrium neutrons and protons and their 
ratios for the same reaction are shown in Fig. [4.6| . When varying the compressibility 
K from 210 MeV (open squares) to 380 MeV (filled circles), it is seen that although 
the yields of both protons and neutrons increase, their ratio remains similar for all 
impact parameters. This is simply because the effects of an isospin-symmetric EOS 
on neutrons and protons are identical. 

The experimental neutron-proton cross section is about three times the neutron- 
neutron (proton-proton) cross section in the energy range studied here. Setting the 
two cross sections equal (fancy squares), it has been found that the proton and neu- 
tron yields and their ratios change by less than 10% even in peripheral collisions of 
^^^Sn +^^^ Sn. This result is also easy to understand since both colliding nucleons 



have the same probability to gain enough energy to become unbound [74|. Thus, the 



in-medium, isospin-dependent nucleon-nucleon cross sections do not affect much the 
ratio i?„/p(i?kin)- In the absence of Coulomb and symmetry potentials their ratio is 
thus almost independent of the nucleon kinetic energy and has a constant value of 
about 2.1 ± 0.3 in both central and peripheral collisions of ^'^^Sn +^^^ Sn. 

Including the Coulomb and the asymmetric term of the EOS in Eq. (^4.2| ), one 
can then study the effects of the symmetry energy, esyni(p), since the Coulomb ef- 
fect is well-known. The symmetry potential has the following effects on preequilib- 
rium nucleons. First, it tends to make more neutrons than protons unbound. One 
therefore expects that a stronger symmetry potential leads to a larger ratio of free 
neutrons to protons. Second, if both neutrons and protons are already free, the sym- 
metry potential makes neutrons more energetic than protons. These effects are clearly 
demonstrated in Fig. where the ratios Rn/p{Ekin) are shown as a function of kinetic 



51 



1 - 

5 -+ 



4 - 

3 r 



I ' ' ' ' I ' ' ' ' I 
b; 

K=380 MeV, 



'^^Sn+^^^Sn, 1 gbg 5 fm 



• K=380 MeV 
□ K=210 MeV 



■ i ■ M ■ 



7 gbg 12 fm 



_L 



_J_ 



20 40 60 80 
Eidn (MeV) 



100 



Figure 4.5: The ratio of preequilibrium neutrons to protons as a function of ki- 
netic energy in central (upper window) and peripheral (lower window) reactions of 
^^^Sn -|-^^^ Sn at a beam energy of 40 MeV/nucleon without the Coulomb and sym- 



metry potentials. Taken from [dl 



energy. These results are obtained by using the three forms of F{u) for typical central 
(left windows) and peripheral (right windows) collisions of ^^^Sn+^^^Sn, ^^^Sn+^^^Sn 
and ^^^Sn -|-^^^ Sn, respectively. The increase of the ratios at lower kinetic energies 
in all cases is due to Coulomb repulsion which shifts protons from lower to higher ki- 
netic energies. On the other hand, the different ratios calculated using different F(m)'s 
reflect clearly the effect mentioned above, i.e., with a stronger symmetry potential 
the ratio of preequilibrium neutrons to protons becomes larger for more neutron rich 
systems. Between the central and peripheral collisions the ratios vary smoothly as a 
function of impact parameter. This is shown in Fig. |4.8| where the energy-dependence 
of the ratio Rn/p{E]an) is shown for different impact parameters for the reaction of 
^^^Sn+^^^ Sn using F-^iu). The impact parameter dependence of the ratio of the total 
yield of neutrons to that of protons shown in Fig. |4.9| clearly demonstrates effects of 
the neutrons skins in more neutron-rich system at larger impact parameters. 

It is interesting to note that effects due to different symmetry potentials are seen 
in different kinetic energy regions for central and peripheral collisions. This can 
be understood from the impact parameter dependence of i?„/p(i?kin), shown in Fig. 



478| , and the reaction dynamics. In central collisions, effects of the symmetry poten- 
tial are most prominent at higher kinetic energies. This is because most of finally 
observed free neutrons and protons are already unbound in the early stage of the 
reaction as a result of violent nucleon-nucleon collisions. The symmetry potential 
thus mainly affects the nucleon energy spectra by shifting more neutrons to higher 
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Figure 4.6: The total numbers of preequilibrium neutrons and protons (lower window) 
and their ratio (upper window) as functions of impact parameter for the reaction of 
^^^Sn +^^^ Sn at a beam energy of 40 MeV/nucleon. 
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Figure 4.7: The ratio of preequilibrium neutrons to protons as a function of kinetic 
energy in central (left windows) and peripheral (right windows) reactions using the 
three forms of F{u). Taken from PT| . 
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Figure 4.8: Impact parameter dependence of the ratio Rn/p{Ei^m) as a function of 
kinetic energy for the reaction of -"^^^Sn +-'^^^ Sn using Fs{u). 
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Figure 4.9: The ratio of preequihbrium neutrons to protons as a function of impact 
parameter using the three forms of F{u). 
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kinetic energies with respect to protons. In peripheral coUisions, however, there are 
fewer nucleon-nucleon colUsions; whether a nucleon can become unbound depends 
strongly on its potential energy. With a stronger symmetry potential more neutrons 
(protons) become unbound (bound) as a result of a stronger symmetry potential, but 
they generally have smaller kinetic energies. Therefore, in peripheral collisions effects 
of the symmetry potential show up chiefly at lower kinetic energies. For the two 
systems with more neutrons the effects of the symmetry potential are so strong that 
in central (peripheral) collisions different forms of F{u) can be clearly distinguished 
from the ratio of preequilibrium neutrons to protons at higher (lower) kinetic ener- 
gies. However, because of energy thresholds in detectors, it is difficult to measure 
low energy nucleons, especially neutrons. Furthermore, the low energy spectrum also 
has appreciable contribution from emissions at the later stage when the system is in 
equilibrium. Therefore, the measurement of the ratio Rn/piE]^m) in central heavy-ion 
collisions for nucleons with energies higher than about 20 MeV is more suitable for 
extracting the EOS of asymmetric nuclear matter. The beam energy range where 
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Figure 4.10: The ratio of preequilibrium neutrons to protons as a function of impact 
parameter for the reaction of ^^^Sn +^'^^ Sn at a beam energy of 100 MeV/nucleon. 

the symmetry potential is relevant for heavy-ion collisions depends on the isospin 
asymmetry of the reaction system and the observables to be studied. To observe 
effects of a weak mean-field potential, such as the symmetry potential, one needs to 
use relatively low beam energies so that the dynamics is not dominated by nucleon- 
nucleon collisions. On the other hand, to study the density dependence and to reach 
a stronger mean-field potential the reactions should be energetic enough to achieve 
sufficient compression. Thus, it is necessary to study theoretically the beam energy 
dependence of the isospin effects on preequilibrium nucleons. Preliminary studies 
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on the neutron/proton ratio in the reaction of ^^^Sn +^^^ Sn indicate that the beam 
energy should not be greater than about 100 MeV. To illustrate this, we show in 
Fig. [4.10| the i?„/p(i?kin) ratio as a function of impact parameter for this reaction at 
a beam energy of 100 MeV. Calculations corresponding to three typical cases of the 
symmetry energy, i.e., Csym = 0, 33 MeV, and 33p/po MeV, have been carried out. 
It is seen that differences among results from these calculations are rather small and 
are all close to the N/Z ratio of the reaction system, i.e., 1.48. 

There seems already experimental evidence in heavy-ion collisions for the exis- 
tence of a strong symmetry potential. For example, in experiments which measure 
preequilibrium nucleons, a number of puzzles have been observed, and they can be 
explained at least qualitatively by effects of the symmetry potential. In the follow- 
ing, several of these puzzles are discussed. In experiments of heavy-ion collisions 
around the Fermi energy ||143| , |144| , |145| , |146|| , Hilscher et al. have found that the neu- 
tron/proton ratio, ((A^/Z)free), of preequilibrium nucleons is consistently higher than 
that of the projectile-target system, {Np + Nt)/{Zp + Zt), and cannot be explained by 
the Coulomb effect. For example, in the reaction of ^^C+^^^Ho at a beam energy of 32 
MeV/nucleon, they have measured the neutron and proton spectra at 14° and energies 
between 70 and 100 MeV, and found that the multiplicity of neutrons is larger than 
that of protons by a factor of 1.4 ± 0.2, 1.7 ± 0.3 and 2.4 ± 0.3 for linear-momentum 
transfers of 52%, 73% and 93%, respectively. Therefore, the neutron to proton ratio is 
much higher than that of the reaction system (A^/Z)cs=1.42 in central collisions cor- 
responding to higher linear-momentum transfer. This result cannot be explained by 
the standard Fermi jet model for preequilibrium nucleon emission ||143|| . On the other 



hand, it is in agreement with the BUU predictions discussed above. From Fig. [4.7| and 
Fig. [4.9| , it is seen that with a large symmetry potential the n/p ratio for energetic 
nucleons can be much larger than that of the colliding system. Quantitative compar- 
isons with the experiments will be very useful for extracting the isospin-dependent 
EOS. Semi-quantitative comparisons between theoretical results and some of the data 
have already been made using a simple symmetry potential energy density of Eq. 
( p.24|) within the buu/lv model |7^, |7^. In Fig. [4.11| the relative neutron excess of 
preequilibrium nucleons 5pre-eq = {N — Z) / {N + Z) is shown as a function of time 
for the reaction of ^^C +^^^ Ho at a beam energy of 35 MeV/nucleon and an impact 
parameter of 1.0 fm. In the early stage of the reaction, the neutron excess is seen 
to fluctuate due to violent collisions, but reaches a constant value after about 100 
fm/c. The dashed line at 0.17 is the neutron excess of the reaction system. With- 
out the symmetry and Coulomb potential the final neutron excess of preequilibrium 
nucleons is about the same as that of the reaction system. In the case of including 
the Coulomb but not the symmetry potential (c = 0) more protons are emitted in 
the preequilibrium stage of the reaction due to the Coulomb repulsion. This effect 
has also been observed in other BUU calculations fTO], |7l], Furthermore, the ratio 



{N/Z)fj.ee increases monotonically as strength of the symmetry energy increases . To 
have a ratio larger than the N/Z of the reaction system one needs to use a strength 
parameter c > 20 MeV corresponding to esyin(O) > 32 MeV. More detailed compar- 
isons with data are required to see what form and strength of the symmetry potential 
are needed to reproduce the data. 
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Figure 4.11: The ratio of preequilibrium neutrons to protons as a function of time 
for the reaction of ^^C +^^^ Ho at a beam energy of 35 MeV/nucleon and an impact 
parameter of 1.0 fm. Sohd circles: c = 46 MeV with Coulomb; open squares: c = 32 
MeV with Coulomb; open circles: c = MeV with Coulomb; solid diamonds: c = 
MeV without Coulomb. Taken from fT^ . 



58 



A more quantitative comparison between the recent data from Ganil and the 
LV model calculation ||75|] has been carried out. In Fig. |4.12| , the experimentally 
measured correlations between neutron multiplicities, M^v, in the forward direction 
(6° < 6* < 20°) and fragment charges, Z/, in the reaction of Pb+Au at i?beamM= 29 
MeV are shown as contours. The contour lines separating the areas in different shades 
(from white to dark gray) have the values 0.6, 1, 2, and 3, respectively. Calculations 
using the LV model without the symmetry potential (c = 0) are shown as open symbols 
and are seen to significantly underpredict the neutron multiplicities, especially in 
central collisions. In calculations with a symmetry energy strength parameter c = 20 
MeV (filled circles) the emission of neutrons is enhanced, so the data are better 
described. 




Figure 4.12: Experimental neutron multiplicity versus fragment charge for peripheral 
reactions of Pb+Au at a beam energy of 29 MeV/nucleon. LV model calculations 
with (without) the symmetry potential are given by filled (open) symbols. Taken 
from |]75| . 

Another interesting experimental observation is the ratio of free neutrons and 
protons from central reactions of ^^S+^^^'^^^Sm — ^ F, F+n,p, where F denotes nuclear 
fragments, at Eheam/A = 26 MeV as shown in Fig. |4.13| . We note that protons emitted 
at velocities higher than the projectile velocity are mainly preequilibrium particles. 
Several interesting observations can be made from these data. First, the emission of 
protons in the neutron-richer system (S +^^^ Sm) is suppressed although both systems 
have the same number of protons to begin with. As we have discussed in detail in the 
previous section, protons feel an attractive symmetry potential, and the emission of 
high energy protons is thus suppressed with respect to neutrons. Correspondingly, the 
emission of high energy neutrons is enhanced, and this is indeed consistent with the 
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neutron data at higher kinetic energies. Second, even for neutrons with lower kinetic 
energies the ratio of free neutrons from the two systems is about 1.23 and is much 
larger than the ratio (1.12) of neutrons of the two systems. More recently, Schroder 
et al. have studied systematically the spectra of preequilibrium neutrons and protons 
in both isospin symmetric and asymmetric collisions ||147|| . It is found that the n/p 
ratio of preequilibrium nucleons does not scale with the N/Z ratio of the combined 
system. In particular, a preliminary analysis of the experimental data indicates that 
the ratio (n/p)48ca+ii2sn/ ('^/p)4oca+ii2sn at a beam energy of 35 MeV/nucleon is about 
4 to 12 depending on the impact parameter. A similar observation has been made 
by Hilscher et al. Both results are thus consistent with the existence of a strong 
symmetry potential. 
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Figure 4.13: The ratio of free neutrons (open circles) and protons (solid diamond) in 
central reactions of ^^S -)-i44,i54 ^ function of the kinetic energy of the emitted 



particle, at a beam energy of 26 MeV/nucleon. Taken from ||144 



Enhanced production of free neutrons than that predicted by the standard sta- 
tistical model have also been reported in Refs. [0, |148|| . In Ref . ||149|| it was pointed 
out that within the experimental and model uncertainties the production of neutron- 
deficient residues is much higher than the statistical model predictions. Because of 
charge conservation, this result is thus equivalent to an enhanced neutron emission 
as shown in Refs. 



This is again consistent with the existence of a strong 
symmetry potential. To extract quantitative information about the nuclear symmetry 
potential requires, however, more detailed studies using the transport model. 

To summarize, we have shown that in the isospin-dependent BUU model for heavy- 
ion collisions at intermediate energies the calculated ratio of preequilibrium neutrons 
to protons is very sensitive to the symmetry potential. Longstanding puzzles observed 
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in preequilibrium neutrons and protons can be understood qualitatively by consider- 
ing effects of the symmetry potential. More quantitative comparisons with data are 
still needed to constrain the form and strength of the symmetry potential. 



4.3 Isospin equilibrium as a probe of nuclear stop- 
ping power 

The degree of stopping in heavy-ion collisions is an important quantity in determining 
the outcome of a reaction. It is usually studied by measuring the rapidity distribution 
or the asymmetry of the nucleon momentum distribution. If we were able to tag the 
nucleons from projectile and target, then this task would be much easier. By using 
nuclei with different N/Z ratios, such a tag can be provided. Recently, it has thus 
been suggested that the degree of equilibration in the isospin degree of freedom also 
provides a means to probe nuclear stopping in heavy ion collisions ||150|| , as illustrated 
in Fig. |4.14| . 



The pictures on the left show the nucleon distributions in coordinate space, while 
that on the right are the rapidity distributions for the three scenarios of pile-up and 
bounce-back, stopping and mixing (isospin equilibrium), and transparency. It is seen 
that for a projectile and a target with very different N/Z ratios, a comparison of the 
rapidity distribution of N/Z before (shown in top picture on the right) and after the 
collision can give direct information about the degree of stopping between the target 
and projectile ||129| , |15CI|| . This method thus allows one to study whether there is 



a transition from full stopping to transparency as the beam energy varies. Also, if 
isospin equilibrium can be reached in the collisions, it is then possible to determine 
its time scale relative to that for thermal equilibrium. Moreover, one can study the 
dependence of isospin equilibrium on the nuclear EOS and the in-medium nucleon- 
nucleon cross sections. We note that this method is restricted to beam energies that 
are below the pion production threshold; so the total nucleon charge is conserved. 

These ideas have been studied in the isospin-dependent qmd model under various 
conditions as shown in Fig. [4.15| . In the left window central collisions of ^°Cr +'^^ Ca 
and ^''Cr Cr at ii^boam/^ = 1-0 GeV are compared, and it shows clearly that there 
is a transparency in the asymmetric system. In the right window two calculations for 
the asymmetric system for different beam energies are shown, and it is seen that even 
at 150 MeV/nucleon the asymmetric system shows transparency. Bass et al. further 
show that the signature for transparency seen in the N/Z ratio is not altered by cluster 
formation. On the other hand, the stopping power is affected by the magnitude of 
the nucleon-nucleon cross section. As expected, increasing the in-medium nucleon- 
nucleon cross section by a factor of 5 results in a transition from transparency to full 
stopping. 

The idea of using isospin to study the nature of stopping in nuclear collisions 
has been introduced earlier in Ref. ||129|| to explain the angular distribution of isotope 
ratios observed in experiments by Yennello et al. ||151| , |152| , |153|| . In these experiments 



the isotopic composition of intermediate mass fragments from central collisions of 
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Figure 4.14: Illustration of using the rapidity distribution of neutron/proton ratio as 



a probe of nuclear stopping power. (A similar illustration can be found in Ref. [|15CI|| . 
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Figure 4.15: The neutron to proton asymmetry versus rapidity predicted by the qmd 



model. Taken from Ref. 150 
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Figure 4.16: Isobaric ratios from central collisions plotted as a function of N/Z ratio 
of the combined target and projectile system at £'beam/^=35 and 45 MeV. Taken 



from Ref. 153 . 
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^OQ, ^°Ar and ^°Ca with ^^Fe and ^^Ni at EbeamM=25, 35, 45 and 53 MeV was 
studied at nscl/msu and tamu. It was found that at i?bcamM=25 and 33 MeV 
the isotopic ratios ^Be/''Be, ^^B/^°B and ^^C/^^C increase hnearly with increasing 
{N/Z)cs ratio of the combined target and projectile system, but are independent of 
the N/Z ratio of the target or projectile. Shown in the left window of Fig. |4.16| are 
typical results of reactions at 33 MeV/nucleon, which indicate that the isospin is 
equilibrated in the composite system formed in these reactions before the emission of 
fragments. 

The most striking and unexpected feature was observed from the isotopic ratios 
in central collisions at E'beamM = 45 and 53 MeV. A typical result at Ebeam/A =45 
MeV is shown in the right window of Fig. 4.16| . It is seen that the isotopic ratios 
depend on the N/Z ratio of the target and projectile in reactions with target-projectile 
combinations having the same {N/Z)cs ratio, such as "^^Ca Fe and ^''Ar Ni. 
In ref. ||153|| , similar results are also shown for other isotope ratios. 



Moreover, data at very forward and backward angles show that the isotope ratios 
do not simply depend on {N/Z)cs- Instead, light fragments at backward angles are 
seen to have a much stronger dependence on target, while at forward angles they 

depend more on (iV/^) projectile- These results demonstrate that the isospin degree of 
freedom in reactions at i?beamM = 45 and 53 MeV is not equilibrated during the 
time when fragments are emitted. Therefore, a transition from isospin equilibration 
to non-equilibration is observed as the beam energy increases from below to above 
the Fermi energy. 

The above observation has profound implications on the reaction mechanism lead- 
ing to multifragmentation. It not only establishes the relative time scale for multi- 
fragmentation in these reactions but also indicates that the assumption of isospin 
equilibrium taken for granted in various statistical models for nuclear multifragmen- 
tation at intermediate energies is not valid. Indeed, a statistical model study was 
made and failed to show any entrance channel effect [ |151|| . Calculations using an 
intranuclear cascade code ISABEL ||155|| show that the N/Z ratio of the residue is very 



close to that of the initial combined system ||151|| and thus also fails to reproduce 
those features observed at i^beam/^ = 45 and 53 MeV. Although the exact origin of 
this failure is not clear, one expects that the reaction dynamics at these relatively 
low energies cannot be described by the nucleon-nucleon cascade alone, and should 
include also nuclear mean-field potential. 

The experimental observation discussed above can be explained by the isospin- 
dependent BUU model outlined in section [4.1|. Calculations based on this model have 



been performed over a range of impact parameters. For peripheral collisions it shows 
a memory of the initial target and projectile, which is, however, gradually lost as the 
collisions become more central. A calculation at b=0 thus gives the most interesting 
test of any non-equilibrium effect. In Fig. |4.17|, we show the energy dependence of 
the density contours in the reaction plane at t= 200 fm/c in head-on collisions of 
Ar+Ni at E'beam/^ = 25, 35, 45 and 55 MeV. The solid contours with p = po/8 
correspond essentially to bound composite systems or heavy residues formed in the 
reactions. The dotted contours with p = 0.05 po are free nucleons emitted mainly 
before thermal equilibrium is reached. The most interesting feature in Fig. [4.17| is the 
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Figure 4.17: Density contours in the reaction plane at t= 200 fm/c in head-on coUi- 



sions of Ar+Ni at ^beam/^ = 25, 35, 45 and 55 MeV. Taken from Ref. pl . 
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formation of two heavy residues in the reaction at i?beam/^ = 55 MeV. For reactions 
at £^beam/^ = 25, 35 and 45 MeV, these heavy residues do not break up even after 
300 fm/c. 
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Figure 4.18: The neutron and proton numbers in the residues on the left {Z < 0) 
and right {Z > 0) side of the origin. Sohd hues are the proton number from the 
projectile, while dot-dashed lines are the neutron number from the projectile which 
moves towards the right. Dashed lines are the proton number from the target, while 
dotted lines are the neutron number from the target which moves towards the left. 



Taken from Ref. [12S 



Studies on similar systems at much lower energies show that the isospin degree 
of freedom is one of the fastest to equilibrate ||156|| . To see if the heavy residues 



observed above is also in isospin equilibrium, we show in Fig. |4.18| the neutron and 
proton numbers in the residues on the left {Z < 0) and right {Z > 0) side of the 
origin in head-on reactions of Ar-|-Ni at i^beam/^ =25, 35, 45 and 55 MeV. The 
solid and dot-dashed lines are, respectively, the proton and neutron numbers from 
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the projectile, which is incident from the left. The dashed and dotted lines are, 
respectively, the proton and neutron numbers from the target, which moves from 
the right. It is seen that the neutron and proton numbers on both sides not only 
decrease but also fluctuate. The decreases is mainly due to nucleon-nucleon collisions 
and particle emissions, while the fluctuation is due to both the restoring force from 
the mean- field potential and nucleon-nucleon collisions. At ii^beam/^ = 25 MeV the 
neutron and proton numbers on the two sides become very close to each other and the 
amplitude of oscillation is rather small by the time of 300 fm/c. This indicates that 
the heavy residue is very close to isospin equilibrium, i.e., the proton and neutron 
distributions are independent of space-time. The damping of the oscillation is faster 
at E'beam/^ = 35 McV so the particle distribution reaches isospin equilibrium also 
sooner. The isotopic composition of fragments emitted from the residues in these low 
energy reactions after about 300 fm/c would therefore essentially refiect the {N/Z)cs 
ratio of the initial composite system, and there is little forward-backward asymmetry. 
These features are in good agreement with those found in the data of i^beam/^ =25 
and 35 MeV P2|, [mHI^l- 



At higher energies, such as -EbeamM = 45 and 55 MeV, there is little oscillation in 
the overlapping region between the target and projectile. This is mainly because the 
incoming momenta of projectile-nucleons and target-nucleons are very large so the 
mean-field potential cannot reverse the directions of motion of many nucleons during 
a relatively short reaction time. As a result, there exists a large isospin asymmetry 
or non-equilibration at these two energies. In particular, on the left side of the origin 
the N/Z ratio of the residue is more affected by that of the target while on the right 
side it is more affected by that of the projectile. However, the N/Z ratios on both 
sides are not simply those of the target and projectile but a combination of the two 
depending on the complicated reaction dynamics. In the case of i^beam/^ = 55 MeV, 
at the time of about 200 fm/c the heavy residue has broken up into two pieces with 
some longitudinal collectivity. The forward moving residue has an average N/Z ratio 
of about 6/5 and an excitation energy of about 8.6 MeV/nucleon, while the backward 
moving residue has an average N/Z ratio of about 16.5/14 and an excitation energy of 
about 6.8 MeV/nucleon. Both residues are thus in approximate thermal equilibrium 



in their own center of mass frame ||129|| but not in thermal equilibrium with each 
other. 

The signature of nuclear transparency from studying the isotopic ratios is consis- 
tent with that from studying the kinematic observables, such as the scaled quadrupole 
moment Qzz/^res defined as 

Q.JAUt) = iMrcs/ ^{2pl-pl-pl)f{r,p,t), (4.3) 

where f{f,p,t) is the Wigner function from the BUU model calculations and Aj-^s is 
the mass of the residue. A zero value of Qzz is a necessary but not a sufficient condi- 
tion for thermal equilibrium, and a positive value of Qzz signals transparency. The 
quadrupole moments of heavy residues formed in the reactions at Eheam/^ = 25 and 
35 MeV are about zero by the time of about 300 fm/c, indicating the establishment 
of thermal equilibrium. Since the heavy residue formed in the reaction at Ebeam/A = 
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55 MeV has already broken up into two pieces at about 200 fm/c, it is interesting to 
examine more closely the evolution of Qzz/^res in this reaction, which is shown in Fig. 
4.19| for particles at local densities larger than po/8 foi' the reactions at E/A= 45 and 
55 MeV. It is seen that Qzz has a small positive value when the system breaks up, 




Figure 4.19: The scaled quadrupole moment of heavy residues formed in head-on 
collisions of Ar+Ni at E'beam/^ = 45 and 55 MeV. Taken from Ref. 



indicating thus a transparency in the collisions. This finding is also in good agreement 
with the kinematic signature of transparency found recently at Ganil for the same 
reaction system at similar beam energies using the indra detector ||157|| . We note 



that in recent calculations using the Antisymmetrized Molecular Dynamics (amd) 
transparency was found at 35 MeV/nucleon for head-on collisions of ^°Ca +^*' Ca 



158|| . Moreover, the signatures of transparency from both the isotopic ratios and the 



quadrupole moments are also in good agreement with earlier BUU predictions on the 



systematics of nuclear stopping power |159 



Effects of the symmetry potential on isospin equilibrium have been studied at lower 
beam energies where isospin equilibrium is usually reached. Within the LV model 
using the simple symmetry potential energy density of Eq. (|]2D @, it has been 



found that a stronger symmetry potential enhances the degree of isospin equilibrium. 



Shown in Table [Ol are the values of N—Z in the quasitarget (QT) and quasiprojectile 
(QP) formed in the reaction of ^^Ca-|-^°Ca at an impact parameter of 7 fm and a beam 
energy of 15 MeV/nucleon. For comparisons, results from reactions of a symmetric 
system ^°Ca Ca are also listed in the table. In both reactions the total mass 
and charge numbers are the same. The difference in the values of N — Z for the 
symmetric system is completely due to the numerical fluctuations of the calculations. 
In the asymmetric reaction the initial value of iV — Z is 16 and for the projectile 
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Table 4.1: Values of N-Z for the quasitarget(QT) and quasiprojectile (QP) 



Reaction/N 


-Z 


c=0 


c= 


=20 


c=28 


s^Ca Ca 


QT 


3.74 


4. 


64 


5.02 




QP 


8.76 


5. 


82 


4.92 


48ca +48 Ca 


QT 


6.52 


5. 


78 


5.52 




QP 


6.96 


5. 


84 


5.12 



and target, respectively. It is seen that significant stopping is achieved even in the 
case of no symmetry potential. However, without the symmetry potential isospin 
equilibrium cannot be reached, and the degree of isospin equilibrium increases with 
the strength of the symmetry potential. 

In summary, isospin equilibration in heavy-ion collisions can be used to study the 
nature of stopping. The degree of isospin equilibrium depends sensitively on both the 
in-medium nucleon-nucleon cross section and the EOS of asymmetric nuclear matter. 
In both experiments and BUU model calculations a transition from isospin equilibrium 
at lower energies to nonequilibrium at higher energies has been observed. The isospin 
nonequilibrium signals the onset of nuclear transparency. For head-on collisions of 
the system Apj-o + Atar ~ 40 + 58 the transparency starts at a beam energy as low as 
45 MeV/nucleon. 



4.4 Isospin dependence of nuclear collective flow 



Nuclear collective flow in heavy-ion collisions at intermediate energies has been a 
subject of intensive theoretical and experimental studies during the last decade; for 
a general introduction and overview see Ref. | |160| |. The study of the dependence of 
collective flow on entrance channel parameters, such as the beam energy, mass number 
and impact parameter, have revealed much interesting physics about the properties 
and origin of collective flow. In particular, from studying the beam energy dependence 
it has been found that the transverse collective flow changes from a negative one to a 
positive one at an energy -Ebah defined as the balance energy, due to the competition 
between the attractive nuclear mean-field potential and the repulsive nucleon-nucleon 
collisions [[161], [162|, [T6|, |T6|, |T65|, |T66|, |T68|, |T6|. The balance energy has been 



found to depend sensitively on the mass number, impact parameter and properties of 
the colliding nuclei, such as the thickness of their surfaces |p.21 



theoretical studies using transport models (for a review see e.g. ||131| , [173| , |134 



Furthermore, detailed 
have 
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shown that both the strength of transverse flow and the balance energy can be used 
to extract information about the nuclear equation of state and in-medium nucleon- 
nucleon cross sections |8|, |8|, |T2|, |T36|, |l7g, [ITT], |I7g, |T74|, |I7|, ^ . 



200 



[] 

J 

■a 100 b--'- I 



70 
50 



T — I I I I I I I I I I I 1 1 1 1 1 1 1 1 iiii|iiii|iiii|iiii| 1 1 — I — r 



J I I I I I I I I I I I I I I I I I I I llllllllllllllllllll I I I L. 



20 30 50 70 100 
A^+Ap 



200 



Figure 4.20: The balance energy as a function of the combined mass. The experi- 
mental data are compared with BUU calculations using a soft EOS and a reduction of 
in-medium nucleon-nucleon cross section by 0% (crosses), 10% (diamonds) and 20% 
(squares) are compared with the data (filled circles). Taken from Ref. | |120| . 



Shown in Fig. [4.20| is the mass dependence of the balance energy measured by 
Westfall et al. in comparison with the predictions of Klakow et al. using the BUU 
model [|121|| . It is seen that BUU calculations reproduce the measured scaling of 



balance energy with mass, but underpredict the values of balance energy. Using an in- 
medium cross section of Eq. ( p.6|) with a parameter of a = —0.2 the BUU calculations 
can better reproduce the data. This results thus indicate that the collective flow and 
balance energy depend strongly on the in-medium nucleon-nucleon cross section. 

Nuclear transverse flow is expected to be isospin-dependent for several reasons. 
First, it is well-known that nucleon-nucleon collisions cause repulsive collective flow, 
and this effect is roughly proportional to the number of collisions in the overlapping 
region. While the number of particles in this region might not be so strongly isospin- 
dependent, the number of collisions during the reaction of two neutron-rich nuclei is 
smaller. This is due to both the smaller neutron-neutron cross sections (see Chap- 
ter ^ and the larger Pauli blocking of neutron-neutron collisions in these reactions. 
This effect is stronger in peripheral collisions of neutron-rich nuclei where two thick 
neutron skins overlap during the reaction. Second, the Coulomb potential also causes 
repulsive scatterings, and this effect is obviously weaker in the neutron-richer system. 
Third, the isospin-independent part of the nuclear mean-field potential is attractive at 
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densities up to about 2po and 4po, respectively, for the stiff {K = 380 MeV) and soft 
{K = 200 MeV) EOS. Since the mean field effect is approximately proportional to the 



total surface area of the colliding system ||121||, collisions of neutron-rich nuclei are ex- 



pected to be influenced by stronger attractive mean-field potentials due to their more 
extended surfaces. Finally, the generally repulsive symmetry potential is stronger 
for collisions of neutron-rich nuclei. However, since the magnitude of the symmetry 
potential is smaller than the isospin-independent mean-field potential, effects of the 
symmetry potential on the collective flow are expected to be small. Although a more 
quantitative study on the relative importance of the above mechanisms remains to be 
worked out, one expects that the overall isospin effect is that the neutron-richer sys- 
tem will have a stronger negative flow below the balance energy and a weaker positive 
flow at higher energies. Moreover, the relative importance of different contributions 
to the collective flow depends on the beam energy. As the beam energy increases 
the repulsive nucleon-nucleon collisions become dominant and effects of the neutron 
skin is thus less important. Also, the isospin dependence of the nucleon-nucleon cross 
sections becomes weaker at higher energies as discussed in Chapter ^. It is therefore 
understandable that the isospin effects on the collective flow decrease as the beam 
energy increases. These discussions thus show that both the strength of flow and the 
balance energy are expected to be isospin-dependent. 

The isospin dependence of collective flow and balance energy has recently been 
studied theoretically with isospin-dependent BUU models |p.30| , |177|| and experimen- 
tally by Pak et al. at nscl/msu [|178|, |179|| . Within the isospin-dependent BUU model 



the standard transverse momentum analysis ||170|] (see also ||160|| ) are carried out for 
two reaction systems, ^'^Cr+^^Ni and ^^Ca-|-^^Fe, which have the same mass number 
of 48 -f 58 but different neutron/proton ratios of 1.04 and 1.30, respectively. Typical 
results for central collisions at an impact parameter of 2 fm and beam energies of 50, 
60 and 70 MeV/nucleon are shown in Fig. [4.21| . At a beam energy of 50 MeV/nucleon, 
the transverse flow in the reaction of ^^Ca Fe is still negative while that in the 
reaction of ^^Cr Ni becomes already positive. The difference disappears at beam 
energies above 70 MeV/nucleon. 

More quantitatively, one can compare the flow parameter F deflned as the slope of 
the transverse momentum distribution at the center of mass rapidity ?/cm- The beam 
energy dependence of the flow parameter for the two reaction systems at impact 
parameters of 2 and 5 fm are shown in Fig. |4.22| . The lines are the least-square 
flt to calculations using linear functions, i.e., F(Ca-l-Fe) = —32.2 + 0.55 i^/A and 
F(Cr + Ni) = -23.9 + 0.48 E/A at 6 = 2 fm; and F(Ca + Fe) = -35.9 + Q.22E/A 
and F(Cr + Ni) = -23.2 + 0.18 E/A at 6 = 5 fm. It is seen that in both central and 
peripheral collisions the neutron-rich system ^'^Ca-l-^^Fe shows systematically smaller 
flow parameters, indicating thus a stronger attractive interaction during the reaction. 
The effect is more appreciable in peripheral collisions - as expected. Consequently, 
the balance energy in '^^Ca Fe reaction is higher than that in the reaction of 
'^'^Cr Ni by about 10 to 30 MeV/nucleon. The difference between flow parameters 
in the two systems decreases as the beam energy increases and flnally disappears as 
the beam energy rises far above the balance energy. 

The experimental investigation of the isospin dependence of transverse collective 
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Figure 4.21: The transverse momentum distribution in the reactions of '^^Cr Ni 
and ^^Ca Fe at an impact parameter of 2 fm and beam energies of 50, 60 and 70 
MeV/nucleon, respectively. Taken from Ref. [ [130|| . 
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Figure 4.22: Flow parameters for the reactions of ^^Cr Ni and ^^Ca Fe at 



impact parameters of 2 and 5 fm. Taken from Ref. [130|. 



flow has been recently carried out by Pak et al. et al. at nscl/msu. Since one 
cannot identify the direction of flow using the transverse momentum analysis in ex- 
periments, the absolute value of flow parameter is usually extracted. Shown in Fig. 
4.23| are the flow parameters of particles with charge Z = 1, Z = 2 and Z = 3 as 



functions of reduced impact parameter h/h-ca_s,y^ for the reactions of ^'^Fe Fe and 
^^Ni Ni at a beam energy of 55 MeV/nucleon. Note that at this beam energy 
flow is still dominated by the attractive mean-field potential and is thus negative - 
if the same convention as in the BUU calculations is used. It is seen that the flow 
parameter for the neutron-richer system is consistently higher and is in agreement 
with the BUU predictions. Pak et al. have also studied the flow parameter as a 
function of the isotope ratio of the composite projectile plus target system for three 
different fragment types from three isotopic entrance channels. Shown in the upper 
window of Fig. |4.24| is the mean transverse momentum in the reaction plane ver- 



sus the reduced cm. rapidity for Z=2 fragments from impact-parameter-inclusive 
^^Mn Fe collisions at 55 MeV/nucleon. The flow parameter extracted for the 
impact-parameter-inclusive events is plotted in the lower window of Fig. |4.24 as a 



function of the ratio of neutrons to protons of the combined system {N/Z)cs- The 
flow parameter increases linearly with the ratio {N/Z)cs for all three types of par- 



ticles. Fig. [4.25| shows the measured energy dependence of the flow parameter for 
midcentral collisions. The curves are drawn to guide the eyes. In the whole energy 
range the flow parameter shows dependence on the isospin asymmetry of the colliding 
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Figure 4.23: Flow parameters for the reactions of ^'^Fe Ar^^ Fe and ^*Ni Ni as 
a function of the reduced impact parameter at a beam energy of 55 MeV/nucleon. 
Taken from Ref. [[L7S 
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Figure 4.24: (Upper window) Mean transverse momentum in the reaction plane versus 
the reduced cm. rapidity for Z=2 fragments from impact-parameter-inclusive 
Fe collisions at 55 MeV/nucleon. (Lower window) Isospin dependence of the flow 
parameter for inclusive collisions at a beam energy of 55 MeV/nucleon. Taken from 



Ref. [ITTa 
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Figure 4.25: Energy dependence of the flow parameter for midcentral collisions of 
58pg _^58 pg g^j^^ 58j^j _^58 j^j^ respectively. Taken from Ref. fTTSI . 
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system. The minima of the curves correspond to the balance energies where the flow 
changes sign. Because of the isospin dependence of nucleon-nucleon cross sections, 
there are less number of nucleon-nucleon collisions in the neutron-rich system. Below 
the balance energy attractive mean-field potentials dominate. There fewer collisions 
mean less kinetic pressure and help generate more attractive flow in the neutron-rich 
system. Above the balance energy, where the repulsive nucleon-nucleon collisions 
dominate, less collisions in the neutron-rich system results in a smaller flow parame- 
ter. For both systems studied in the experiment the balance energy occurs at about 
-E'beamM = 90 McV. Thcse features in the measured energy dependence are also in 
quantitative agreement with those predicted in Fig. |4.22| using the BUU model. 

The comparison of theoretical predictions on collective flow with the experimental 
data is shown in Fig. |4.26| . In the upper window the measured and predicted bal- 
ance energies are compared for both systems, and in the lower window the measured 
and predicted differences in balance energies of the two systems are compared. One 
observes that the difference of balance energies, i.e., the 5i?bai, is well reproduced. 
The model, however, underpredicts the balance energy in central collisions and over- 
predicts the data in peripheral collisions. These discrepancies can be understood as 
follows \TK, TT2, 18(J|. We mentioned earlier in this section that in central collisions 



nucleon-nucleon collisions are important, and a reduction of nucleon-nucleon cross 
sections by about 20% was needed to reproduce the measured balance energy - see 
Fig. ^.20| . Once one includes the isospin dependence, we expect a different reduction 
of the in-medium nucleon-nucleon cross sections is needed to reproduce the data. This 
remains to be studied. In the case of peripheral collisions, the mean-field potential 
has a more important effect on flow than nucleon-nucleon collisions. In particular, 
the repulsive momentum- dependent mean- field potential [9^, ^ |9^, which is not 



included in the calculations presented here, is expected to reduce the balance energy 
in heavy ion collisions at large impact parameters | |180 |. 

It is encouraging that both model calculations and experiments show that the 
transverse flow is sensitive to the isospin asymmetry of the colliding systems. Studies 
using the isospin-dependent BUU model indicate, however, that effects of the symme- 
try potential on the transverse flow and balance energy are not appreciable, because 
they are small compared to those of the isospin- independent part of the nuclear EOS. 
By neglecting the Coulomb interaction, the balance energies for both systems in- 
crease, but the difference remains about the same in central collisions and becomes 
larger in peripheral collisions. Thus, the observed isospin dependence of the transverse 
flow and balance energy is mainly due to the isospin-dependent in-medium nucleon- 
nucleon cross sections. The isospin-dependent BUU model reproduces the measured 
isospin dependence of the balance energy, but fails to reproduce the magnitude of the 
balance energy. It thus indicates the need of a reduced in-medium nucleon-nucleon 
cross section and a momentum-dependent mean-field potential. More detailed work 
in refining the comparison between model calculations and the data remains to be 
done to infer useful information about the in-medium nucleon-nucleon cross sections. 
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Figure 4.26: Upper window: comparison of the measured and calculated balance 
energies for collisions of ^^Fe Fe and ^^Ni Ni as functions of impact parameter. 
Lower window: comparison of the difference in balance energies of the two systems 
as a function of impact parameter. Taken from Ref. [p.79||. 
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4.5 Isospin dependence of total reaction cross sec- 
tions and radii of neutron-rich nuclei 



We discuss, in this section, the isospin dependence of total reaction cross sections 
of neutron-rich nuclei. The total reaction cross section provides direct information 
about the size of nuclei involved in the reaction, and has been one of the standard 
techniques for studying the size and structure of radioactive nuclei 0. Theoretical 
analyses of the measured reaction cross sections have been done almost exclusively 
by using the Glauber model or eikonal approximation with various modifications. 
Shown in Fig. [4.27| is the total reaction cross section for ^^C+^^C reactions at beam 
energies from about 10 to 2000 MeV/nucleon |p.81|| . The model works surprisingly 
well in the whole energy range, although it should not be an appropriate model at 
low energies. Indeed, recent Glauber model calculations underpredict significantly 
the interaction cross sections for ^^Li + C reactions at beam energies below about 100 
MeV/nucleon ||182|, |183||. Furthermore, one expects that both the nuclear EOS and the 




Figure 4.27: The energy dependence of the total reaction cross section for the reaction 
of C. Taken from Ref. pl|] . 

in-medium nucleon-nucleon cross sections, which one has learned from intermediate 
energy heavy-ion collisions, should have effects on the total reaction cross section. In 
the following, we shall first review the observed isospin dependence of total reaction 
cross sections, and then comment on the feasibility of applying isospin-dependent 
transport models to study the total reaction cross section. As we shall show that 
results from these studies have already demonstrated the possibility of learning the 
nuclear EOS, in-medium nucleon-nucleon cross section, and the size of nuclei from the 
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total reaction cross sections. 

Several empirical formulae for describing the systematics of the total reaction 
cross sections have been developed based on the Glauber model [|181| , |185| , |184|| . For 
example, Kox's parameterization for the total reaction cross section can be written 

as iig 

aR{E) = nrlfiE), (4.4) 

where f{E) is an empirical expression describing the degree of transparency, and 
depends on the masses of the colliding nuclei and the beam energy. The strong- 
absorption radius parameter tq contains information about the size of both the target 
and projectile. A constant value of tq ~ 1.05 fm independent of the beam energy 
and projectile-target combination has been obtained by analyzing a large amount 
of (Jr{E) data for collisions between /3-stable nuclei at energies from 33 to 2100 
MeV/nucleon. Using the same parameterization for reactions involving neutron-rich 
nuclei Tq has been found to increase with increasing neutron excess in the reaction 
system ||, |8|, |186|] . Shown in Fig. [4.28| are the extracted values of function 
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Figure 4.28: Square of the radius parameter of the interacting nuclei Tq as a function 
of — Z of the reaction system. The open circles, filled diamonds, open diamonds, 
crosses and filled circles are for ^^-^^N +^8 Si, ^e-ss^a +28 gi, ^'^Q +62-68 2n and 
^^O +^8 Si, respectively. Taken from Ref. 
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18^. 



of the neutron excess N — Z oi the reaction system | ]186| . The linear dependence 
of Tq on N — Z can be explained using an extended Glauber model where nucleons 
are assumed to have a Fermi-type distribution with the following isospin-dependent 

PI 

(4.5) 



surface diffuseness for neutrons 



t„ = 2.4 + 10 ■ (5 - 5^) (fm). 



80 



where Sf^ = {{N — Z)/A)g is the relative neutron excess along the P stability line. 




Figure 4.29: Interaction cross section for Na isotopes on carbon at a beam energy of 
950 MeV/nucleon. Taken from Ref. ||187|. 



More recently, Suzuki et al. have measured the interaction cross section of Na 
isotopes (20-23,25-32]\j-g^-j ^ carbou target at a beam energy of 950 MeV/nucleon 
187|| . In Fig. [4.29| the measured cross section (filled circles) is plotted as a function 



of mass number of the Na isotopes. The dotted line is obtained from the expression 
ai = T:(Ri{C)+roA"^)\ (4.6) 



where Ri{C) = 2.61 fm is the interaction radius of ^^C and tq is chosen to reproduce 
the ^^Na C data. It is seen that as the projectile becomes more neutron-rich 
this simple scaling starts to deviate from the data. This result is consistent with 



that observed in Refs. ||184| , |186|| , although different formulae have been used there to 
describe the reaction cross section. The results from the Glauber model are also shown 
in the figure. The solid line is obtained by assuming that the density distribution for 
neutrons is the same as that for protons, except a normalization factor (N/Z). With 
this assumption the model underpredicts the cross section for neutron-rich nuclei 
but overpredicts that for the neutron-poor nuclei. Using the neutron and proton 
density distributions predicted by the RMF theory, the data can be well reproduced 
as shown by the dashed line. The increasing difference between the dashed and solid 
lines reflects the effects due to the larger neutron-skin of Na isotopes. In fact, a 
gradual increase of the neutron-skin up to 0.4 fm for neutron-rich Na isotopes has 
been deduced from fitting the total reaction cross sections using the Glauber model 
T87| . 
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In the optical limit, the Glauber model predicts a reaction cross section 



aR = 27f [1 - T{b)]bdb, 
Jo 



(4.7) 



where T{b) is the transmission function at an impact parameter b. In the standard 
Glauber model, T{b) is evaluated from a convolution of the nucleon-nucleon cross 
section an, with the density distributions of target (T) and projectile (P) nucleons 
in the overlapping region [|185|| , i.e.. 



T(6) = exp (^-Y^^ki J pTk{s)ppi{b- s)ds 



where the indices k, I are used to distinguish neutrons and protons and is the 
nucleon density distribution along the beam direction, i.e.. 



p(V + z'^)dz. 



(4.9) 
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Figure 4.30: The impact parameter dependence of the reaction probability 1 — T{b) 
for the reactions of ^^C C (left window) and ^^Li C (right window) at four 
beam energies. Circles, squares, diamonds, and crosses represent the results for E/A 



43, 75, 400, and 800 MeV, respectively. Taken from Ref. [pi 



In transport models, the transmission function T{b) can be calculated directly by 
simply counting events in which no nucleon-nucleon collision have occurred. This 
approach, however, is hampered by large statistical fluctuations. It is thus numer- 
ically more reliable to determine T{b) by using the observation that the number of 
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nucleon-nucleon collisions encountered by a nucleon in heavy ion collisions at all im- 
pact parameters and for both low and high beam energies can be very well described 
by a Poisson distribution [p.59| , |188| , |189| , |190|| . Therefore, the transmission function 
at a given impact parameter b can be simply expressed by the average number of 
collisions per nucleon N{b), i.e., 



T{b) = exp (-iV(6) 



(4.10) 



The isospin dependence of the reaction dynamics is contained in N{b). However, the 
isospin-dependent transport models have not been used to calculate the total reaction 
cross section. Instead, the normal isospin-independent BUU model was used by Shen 
et al. in calculating the total cross sections [|191|| . In their study neutrons and protons 
are distributed randomly in a sharp sphere with a radius of r = r^A^^^. By varying 
the nuclear EOS, in-medium nucleon-nucleon cross section and the radius parameter rg 
Shen et al. have already demonstrated the usefulness of transport models in studying 
the total reaction cross section. Fig. 4.30 shows the reaction probability 1 — T{b) as a 
function of impact parameter for the reaction of ^^C C and ^^Li C at the four 
beam energies from 43 to 800 MeV/nucleon. It is seen that in central collisions these 
nuclei behave as black spheres. The observed energy dependence of the total reaction 
cross section is due to the energy dependence of partial transparency in peripheral 
collisions. 




Figure 4.31: The energy dependence of the total reaction cross section for three 
different reactions. The data are shown by filled circles while the solid (dashed) lines 
are calculations with the stiff (soft) nuclear EOS. Taken from Ref. ||191|| . 



The calculated and measured energy dependence of the total reaction cross section 



are compared in Fig. |4.31|. The filled circles are the experimental data, while the 
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solid (dashed) lines are the BUU calculations with a stiff (soft) nuclear EOS. Both 
calculations can reproduce the data at beam energies above about 100 MeV/nucleon. 
The dependence of the total cross section to the nuclear EOS appears at lower beam 
energies, but this is exactly where the Glauber-type models underpredict the total 
reaction cross sections ||182| , |183| . 

BUU calculations generally overpredict the measured cross sections at beam en- 
ergies below 100 MeV/nucleon. This might be related to the in-medium nucleon- 



nucleon cross sections used in the model. In Ref. ||191|| a parameterization which fits 
the proton-proton data and has a cutoff of 55 mb at lower beam energies has been 
used for all nucleon-nucleon collisions. If the nucleon-nucleon cross section is mod- 
ified in medium, then the total reaction cross section is expected to be affected as 
well. This effect is shown in Fig. [4.32| , where it is seen that for all incident energies 
a 20% variation of the nucleon-nucleon cross section results in about 15% variation 
of the total reaction cross section. This variation appears to be larger than that 



found in Ref. |192] using the Glauber model where a 7% reduction has been found 



by reducing the effective, energy-independent nucleon-nucleon cross section 40 mb by 
20%. Furthermore, a 20% reduction the BUU model can better reproduce the data. 
This seems to be consistent with that from studying the balance energy [OTi] and 
the impact parameter dependence of the flow parameter ||122|| . However, one should 



note that the 20% reduction of the proton-proton cross section used in Ref. ||191 
for all nucleon-nucleon collisions corresponds to a larger reduction of the in-medium 
nucleon-nucleon cross section if its isospin dependence is taken into account. 

Transport models may also be useful for studying the isospin-dependence of the 
radius parameter rg. Keeping the rg parameter fixed at 1.33 fm for ^^C the calculated 
total reaction cross section for ^^Li +^'^ C has been found to be very sensitive to the 
tq parameter for ^^Li. Shown in Fig. [4.33| are the tq dependence of the total reaction 
cross section. As one expects, the total reaction cross section increases significantly 
with increasing for ^^Li. A radius parameter tq = 1.60 fm is needed to fit the cross 
section data. 

From the above discussion it is clear that transport models are useful for studying 
the total reaction cross section, especially at low beam energies where Glauber-type 
models are not valid. So far only the isospin-independent BUU model has been used 
to study the total reaction cross section and its dependence on the nuclear EOS, in- 
medium nucleon-nucleon cross section, and the size of neutron-rich nuclei. It will be 
interesting to repeat such studies using the isospin-dependent transport models. 



4.6 Isospin dependence of subthreshold pion pro- 
duction 

The multiplicity ratio /n^ of charged pions in heavy-ion collisions is closely related 
to the ratio N/Z in the participant region. A systematic study of the relation be- 
tween the two ratios are useful not only for extracting information about the nucleon 
distributions but also for studying the production mechanism of subthreshold pions. 
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Figure 4.32: Effects of in-medium nucleon-nucleon cross sections on the energy de- 
pendence of total reaction cross section for the reaction of ^^C C. The data are 
plotted with filled circles. All calculations are carried out with the stiff nuclear EOS. 
The solid line is obtained using the free-space pp cross section; while the dashed and 
dotted lines are calculated with a reduction of the cross section by 20% and 40%, 



respectively. Taken from Ref. ||191 
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Figure 4.33: Effects of radius parameter ro on tlie energy dependence of total reaction 
cross section for the reaction of ^^Li+^^C. The data are plotted with filled circles. The 
calculations are represented by the open symbols and connected by the dashed, solid, 
and dotted lines, representing the choices of tq = 1.8, 1.6, and 1.33 fm, respectively. 
All calculations are carried out with the stiff nuclear EOS. Taken from Ref. 



Assuming that pions are produced through the intermediate A(1232) resonance, 
i.e., A^A^ A^A(1232) and A ^ A^ + vr, the ratio n-/n+ can be estimated [|19|, [194 
and is given by 



TT 



+ 5Z'^ + NZ \Z 



This simple scaling formula is supported by data from relativistic heavy-ion exper- 
iments at the BEVALAC ||193|, |195||. Fig. ^4.34| shows the tt~ /-k^ ratio as a function 



of N/Z ratio of the participants in the reaction of Ar+KCl at a beam energies of 
2.1 GeV/nucleon. It is seen that the ti~ /n^ ratio scales well with {N/Z)"^ of the 
participants. At subthreshold energies, however, it has been found that the ratio 
-n^ j-n^ scales approximately linearly with N jZ from in reactions induced by ^^C on 
^Li, 12c, ii'^Sn and i24Sn at ^beamM = 60 - 85 MeV ll9|, For example, at 



-^beam/^ =85 McV the relative ratio for reactions with ^^^Sn and ^^^Sn targets is 



(y-K- 



Sn)/ — (ii^Sn)^1.2±0.3, (4.12) 

which is compatible with the ratio of neutrons in the two targets. Such a linear 
dependence of 7r~/7r+ on NjZ at subthreshold energies is in contrast to the relation 
found at relativistic energies, although the available data is rather limited at both 
energies. It is also interesting to note that an anomalously large ratio of ii~ /tt^ 
7 has been observed in the reaction of ^^C -F^"^ Pb at ^beamM = §5 MeV il97 |. 
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Figure 4.34: The ratio of tt /7r+ as a function of N/Z ratio of the participants in the 
reaction of Ar+KCl. Taken from Ref. |195|. 



This ratio cannot be explained by the neutron-excess in the 208pb target. Thus, the 
relation between tt' /n^ and N/Z of the participants has not been well established 
experimentally at subthreshold energies. The simple scaling law of Eq. ( 4.11 ) thus 



needs to be tested in a wide range of beam energies and N/Z ratios. High intensity 
beams of neutron-rich or radioactive nuclei are useful for this purpose. 

The possibility of using the 7r~/7r"'" ratio to extract information about the distribu- 
tion of neutrons and protons in neutron-rich nuclei has been studied using the Glauber 
model 1 198, 199| , 200 1. Under the assumption that the nucleon density distribution is 



a Gaussian function, i.e., 



2 

p(r)=p(0)exp(-^), (4.13) 



the integration in Eq. (|4.8| ) can be carried out analytically to yield the following 
result 



T{h) = exp(-iV(6)), (4.14) 

with the average number of nucleon-nucleon collisions per nucleon in the reaction of 
A + B at an impact parameter b 

m = '--"p-^^^p-^f<-B . L-l—\ . (4.15) 

In the above a is the isospin-averaged nucleon-nucleon cross section. Similarly, the 
average number of proton-proton, neutron-neutron and proton-neutron collisions can 



87 



be obtained. If one assumes that pions are produced through A(1232) resonances, 
the inclusive 7r+ and vr" cross sections can be written as ||199| , p01|| 



da- 



NA\ 



'ZaZb- 



2tc I bdb ■ exp 

10 



b^ 



a 



Za 



a 



Zb. 



P. 



(4.16) 



dn 



where 



Ps = exp 



2 7r2pjv^(0)pjvs(0)a^^a^^ 
\jna\ ^^a^^b — 



27r J bdb ■ exp ^— 



b' ' 
O-Na + ^Nb . 



■P. 



a(AB-l)p^(0)pB(0)a>| 



exp 



(4.17) 



(4.18) 



is the elastic survival probability, pn- and pz^ are the neutron and proton densities of 
the nucleus i = A or B, and Jna is the amplitude for the reaction N + N ^ N + A. 
Knowing the density distribution of neutrons and protons in both the target and 
projectile, one can then calculate the ratio of -n^ j-n^ from the above relations. It has 
been found that this ratio is rather sensitive to the difference between the neutron 
and proton density distributions. Shown in Table |4.2| are the ratio E{a^ defined as 



E 



(J- 



a- 



(4.19) 



In this calculation ||199|| density distributions from the binding energy adjusted shell 
model [113 are used. Results using two isospin-averaged nucleon-nucleon cross sec- 
tions of 40 mb and 25 mb are compared. It is seen that the ratio E is not so sensitive 
to the cross section a, but is sensitive to the difference between the density distri- 
butions of neutrons and protons. This prediction has attracted some experimental 



interest ||202|| and can be tested in future experiments with high intensity radioactive 
beams. 
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Table 4.2: Normalized difference between tt and tt"*" production cross sections for 
reactions of Li isotopes on ^^C 

Ratio E(f7) ^Li+^^C ^Li+^^C ^Li+^^C ^^Li+^^C 

CT = 40mb 0.1153 0.2221 0.2955 0.3951 
CT = 25mb 0.1143 0.2210 0.2939 0.3927 
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Chapter 5 

Overview and perspectives 



Calculations based on microscopic many-body theories have predicted many interest- 
ing new physics phenomena in isospin- asymmetric nuclear matter. These theoretical 
studies have shown that the equation of state of isospin-asymmetric nuclear matter 
depends quadratically on the neutron excess 5. The coefficient of the quadratic term 
in the EOS of asymmetric nuclear matter, i.e., the symmetry energy, and especially 
its density dependence vary widely among theoretical predictions. In particular, the 
symmetry energy at high densities suffers the most severe uncertainties among all 
properties of dense nuclear matter. 

Theoretical studies also indicate that the isospin-dependent EOS has significant 
influence on the properties of neutron stars as well as on those of radioactive nuclei. 
The softening of the equation of state in isospin-asymmetric nuclear matter is im- 
portant for the prompt explosion of a type ll supernova. The proton concentration 
in neutron stars at j3 equilibrium is almost completely determined by the symme- 
try energy. Consequently, the cooling rate and neutrino emissions of neutron stars 
are affected significantly by the symmetry energy, especially its density dependence. 
The binding energy, radii and deformation of radioactive nuclei are also shown to 
depend critically on the isospin-dependence of the nuclear EOS. The nucleon-nucleon 
cross sections in free space are known to be also strongly isospin-dependent. Nuclear 
many-body theories have predicted that the isospin-dependence of nucleon-nucleon 
cross sections might be significantly altered by the nuclear medium. 

Because of the additional degree of freedom, isospin, new phenomena are predicted 
to occur in isospin-asymmetric nuclear matter. In particular, the chemical instability 
is predicted to occur in a larger configuration space than the mechanical instability. 
Consequently, the spinodal decomposition in asymmetric nuclear matter might be 
triggered by fiuctuations in the isospin asymmetry rather than baryon density. This 
may have significant effects on the mechanism for nuclear multifragmentations in 
nuclear reactions involving isospin asymmetric nuclear matter. The liquid-gas phase 
transition in asymmetric nuclear matter has been predicted to be second-order rather 
than first-order as in symmetric nuclear matter. Moreover, the critical temperature of 
liquid-gas phase transition gradually decreases to zero from symmetric nuclear matter 
to pure neutron matter. 

The predicted new physics phenomena in isospin asymmetric nuclear matter not 
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only motivates further theoretical studies but also calls for dedicated experiments to 
test them. Energetic collisions of stable and radioactive neutron-rich nuclei can cre- 
ate transiently nuclear matter with appreciable isospin asymmetry and matter com- 
pression in the overlap region, thus providing a testing ground for studying isospin 
physics. In fact, many interesting isospin effects have already been observed in heavy- 
ion collisions. For example, isospin non-equilibrium has been observed in the angular 
distribution of isotope ratios. Interesting isospin-dependencc of nuclear multifrag- 
mentation has been observed from comparing reactions of isospin-symmetric systems 
and highly isospin-asymmetric ones. Abnormally preferential emission of preequilib- 
rium neutrons have been observed in many experiments. Very recently, the isospin- 
dependence of nuclear collective flow and balance energy has been first predicted by 
using isospin-dependent transport models and later confirmed by experiments. The 
theoretical understanding of the observed isospin-dependent phenomena has been lim- 
ited to the qualitative level in most cases. For example, isospin-dependent transport 
models, though very useful for predicting and explaining many observations qualita- 
tively, are still unable to explain the data quantitatively. Nonetheless, these models 
are very useful in establishing connections between the experimental observables and 
the isospin-dependent equation of state and in-medium cross sections. For instance, 
it has been shown that the neutron/proton ratio of preequilibrium nucleons is very 
sensitive to the density dependence of the symmetry energy, but not sensitive to the 
in-medium nucleon-nucleon cross sections. On the other hand, the isospin-dependence 
of the collective fiow and balance energy is mainly caused by the isospin-dependence 
of the in-medium nucleon-nucleon cross sections. 

The main goal of studying isospin physics in heavy-ion collisions is to explore the 
properties of nuclear matter in the transition region between symmetric nuclear mat- 
ter and neutron matter and its equation of state. Dedicated radioactive beam facilities 
will provide excellent opportunities to achieve this goal. We anticipate that there will 
be substantial experimental activities in studying isospin physics in the coming years. 
To extract from these experiments more quantitative information about the EOS of 
asymmetric nuclear matter and isospin-dependence of in-medium nucleon-nucleon 
cross sections, it is necessary to combine theoretical and computational techniques 
in the area of nuclear structure, reaction and astrophysics. With a continued dialog 
between theories and experiments, we believe that this goal will be achieved. 
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